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OF THE 


ROYAL ASTRONOMICAL SOCIETY 
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MEETING OF 1959 JANUARY g 
Dr W. H. Steavenson, President, in the Chair 


‘The Pyesident announced the death of William Otto Brunner, an Associate 
of the Society, and paid a brief tribute to his memory, the Fellows standing. 

The President announced that the Council had awarded the Gold Medal 
of the Society to Dr Raymond Arthur Lyttleton for his outstanding contri- 
butions to many branches of theoretical astronomy, in particiilar to the problems 
of the stability of rotating fluid masses, of the constitution of the stars and stellar 
evolution, and of the origin and peculiarities of the solar system. 

The President announced that the Council had awarded the Eddington 


Meda to Dr James Stanley Hey for his outstanding discoveries concerning 
radio emission from the active Sun, radio reflections from meteors, and discrete 
sources of cosmic radio emission. 


‘The election by the Council of the following Fellows was duly confirmed :— 


Dennis lan Gough, University College of Rhodesia and Nyasaland, Private 
sag 167H, Salisbury, 5. Rhodesia (proposed by A. L. Hales); 

I:dward Irving, Department of Geophysics, Australian National University, 
Canberra, Australia (proposed by 5. C. B. Gascoigne); 

Colin Stewart Lindsay Keay, 24 Cashel Street, Christchurch, New Zealand 
(proposed by R. F. Joyce); 

David Layzer, Harvard College Observatory, Cambridge, Mass., U.S.A. 
(proposed by ‘Tl’. Gold); 

Rex Wilheim Mallett, 40 Queen’s Avenue, Warrington, Lancs. (proposed 
by A. C. Sanderson); 

\rthur Jack Meadows, University Observatory, Oxford (proposed by 
H. H. Plaskett); 

“olin Robert Munford, 13 Ladywood Road, Ipswich, Suffolk (proposed by 
R. G. Andrews); 

Alan Eben MacKenzie Nairn, Physics Department, King’s College, New- 
castle-upon-l'yne, 1 (proposed by S. K. Runcorn); and 

*Francis Charles Wykes, 76 Culverden Park Road, Tunbridge Wells, Kent 
(proposed by C. Hunt). 


* Transferring from Junior Membership. 
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MEETING OF 1959 JANUARY 9 
Dr W. H. Steavenson, President, in the Chair 


‘The Pyesident announced the death of William Otto Brunner, an Associate 
of the Society, and paid a brief tribute to his memory, the Fellows standing. 

The President announced that the Council had awarded the Gold Medal 
of the Society to Dr Raymond Arthur Lyttleton for his outstanding contri- 
butions to many branches of theoretical astronomy, in particular to the problems 
of the stability of rotating fluid masses, of the constitution of the stars and stellar 
evolution, and of the origin and peculiarities of the solar system. 

The President announced that the Council had awarded the Eddington 


Med# to Dr James Stanley Hey for his outstanding discoveries concerning 
radio emission from the active Sun, radio reflections from meteors, and discrete 
sources of cosmic radio emission. 


The election by the Council of the following Fellows was duly confirmed :— 


Dennis Ian Gough, University College of Rhodesia and Nyasaland, Private 
Bag 167H, Salisbury, S. Rhodesia (proposed by A. L. Hales); 
Edward Irving, Department of Geophysics, Australian National University, 
Canberra, Australia (proposed by S. C. B. Gascoigne); 
Colin Stewart Lindsay Keay, 24 Cashel Street, Christchurch, New Zealand 
(proposed by R. F. Joyce); 
David Layzer, Harvard College Observatory, Cambridge, Mass., U.S.A. 
(proposed by 'T. Gold); 
Rex Wilheim Mallett, 40 Queen’s Avenue, Warrington, Lancs. (proposed 
by A. C. Sanderson); 
Arthur Jack Meadows, University Observatory, Oxford (proposed by 
»H. H. Plaskett); 
olin Robert Munford, 13 Ladywood Road, Ipswich, Suffolk (proposed by 
R. G. Andrews); 
Alan Eben MacKenzie Nairn, Physics Department, King’s College, New- 
castle-upon-Tyne, 1 (proposed by S. K. Runcorn); and 
*Francis Charles Wykes, 76 Culverden Park Road, Tunbridge Wells, Kent 
(proposed by C. Hunt). 


* Transferring from Junior Membership. 
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The re-election by the Council of the following Fellow was duly confirmed :— 


Arthur Francis O’Donel Alexander, 1 Athelstan Road, Dorchester, Dorset 
(proposed by J. R. Bazin). 


The election by the Council of the following Junior Member was duly 
confirmed :— 


Jeremy Benson Tatum, 1g Montpelier Road, London, N.3 (proposed by 
C. W. Allen). 

Eighty-four presents were announced as having been received since the last 

meeting, including :- 

A. Argoli: Ephemerides ... 1641-1670 (presented by K. S. Jewson); 

J. Gadbury: Ephemerides ... 1719-1740 (presented by K. 5S. Jewson); 
and 

B. Lehnert, (Ed.): Electromagnetic phenomena in cosmical physics. Sym- 
posium No. 6 of the 1.A.U. (presented by the 1.A.U.). 





ON THE SHAPE OF COMETARY TAILS 
B. L. Meek 
(Communicated by H. Bondi) 
(Received 1958 October 24) 


Summary 


On the basis of a particle swarm structure for a comet, the subsequent 
space distribution of particles of different sizes produced within it near 
perihelion is calculated for a series of later instants. ‘The separation of some 
particles from the main part of the comet is produced by radiation pressure (or 
dynamically equivalent effects), and can be represented as a reduction in solar 
attraction, including of course complete reversal. Some properties of the 
resulting locus of particles (of different sizes and effective gravitation) leaving 
the comet at a single instant are discussed; some of the results obtained are 
applied to the particular case of Comet Arend—Roland. The theory provides 
a straightforward explanation of the apparent sunward ‘ spike’ observed for 
this comet in terms of a fan-shaped distribution of emitted particles observable 
only when viewed from positions in or very near its general plane. 





1. Introduction.—The accretion theory of the origin of comets (1) suggests 
that the head of a comet consists of a diffuse swarm of particles which may not 
necessarily be of equal size. In the case of a long-period comet, mutual gravita- 
tion can prevent complete dispersal near aphelion, but near perihelion the Sun’s 
differential gravitational field is so strong that the mutual attractions are negligible 
and the particles move virtually independently. The theory also requires that as 
the comet traverses the small section of its orbit on the perihelion side all the 
particles composing the head pass through the plane of general motion and thus 
tend to bunch more closely together. Because of this, collisions take place within 
the head, producing finer dust and possibly also gas which may, in part at least, go 
to form the tail. We shall consider the paths of these small dust particles after 
the collisions to obtain some indication of their distribution at any later time. 

It is convenient to make some simplifying assumptions. We shall regard the 
motions as taking place in the plane of the comet’s orbit, and initially assume that 
all the fine dust particles are produced simultaneously and instantaneously at a 
point of the comet’s path which we shall take to be at an angular distance « from 
perihelion. Since the comet will be close to perihelion during the interval in 
which we are interested, each particle can be considered as moving independently, 
under the influence of the Sun’s gravitational field and radiation pressure (or 
equivalent effects) only. 

Both of these forces obey the inverse square law. However, radiation pressure 
varies as the surface area of the particle affected, while gravitation varies as its mass, 
that is (roughly) as its volume. Thus radiation pressure will have a greater com- 
parative effect on small particles than large ones, and in the case of a comet some 
particles produced by internal collisions will be so tiny compared with those 
forming the head that they will be swept away from the main body,to form the tail. 


1* 
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If we treat the comet as a single moving point, as is usual when calculating an 
orbit, we can denote its corresponding central acceleration by y/r?. We shall take 
the orbit of this point to be parabolic. It will remain parabolic after the instant 
of collision if we assume that internal collisions have a negligible effect on the 
motion of the comet. At the instant of collision the dust particles will have the 
same speed and direction of motion as the head, and we shall assume these to re- 
main unaltered by the collision* ; however, after this instant the central accelera- 
tion for each particle will be «u/r?, where « is now essentially less than unity. The 
parameter « will be different for different particles, but we shall assume that for 
any one particle it will remain constant in time. ‘The motion of the comet itself 
will be given by the limiting case « = 1; x =o will correspond to a particle for which 
the gravitational force and the radiation pressure exactly cancel; 1>x>0 
represents particles for which gravitational force is the stronger; and o> x those 
for which radiation pressure is the stronger and the net effect an actual repulsion 
from the Sun. In any particular case the visible tail will correspond to values 
of « between unity and some smaller (probably negative) value dependent on 
the density and composition of the tail material. 

‘The problem may now be stated thus: if we consider a set of particles repre- 
sented by values of « in the continuous range 1 > « > — 0, and given the specified 
conditions, what is the locus of these particles at some given time-interval after the 
comet passes the collision point? 

2. Orbits of particles for different values of «.—If the comet (considered as a 
particle) moves in a parabolic orbit with perihelion distance q, the equation of its 
path will be 2g/r=1+cos@ in polar coordinates, where r is the distance of the 
particle from the Sun and @ is measured from the initial line joining the Sun to 
perihelion. 

The orbits of all the dust particles thus pass through the point r = 2qg/(1 + cos), 
#2. We can take the origin of the time scale at the instant the comet passes 
perihelion, and take the instant of collision to be t=t’. We have assumed that at 
t=t' all dust particles have the same velocity as the comet, so the angular momen- 
tum per unit mass (in standard Keplerian orbit notation this is 4 = r?6), which is 
constant for an individual particle, must have the same value for every particle ; 
that is, A is independent of x. Using this fact and the starting conditions at = 
we obtain the orbit of the particle with central acceleration «y/r? in the form 

2q/r =« + cos 6+ (1—x)cos (8—«). (1) 
‘Yo calculate the time ¢ at which a particle starting at 0= at t=?’ reaches a 
point (r, 4) satisfying (1), we have, since h = r°6, 
rt 8 
| hat= | r dB. , (2) 
t’ a 
We now introduce a new variable u, defined by 
u? = 1 — 2g/r(1 + cos #) (3) 
noting that this quantity is dimensionless and identically zero for the case « =1. 
It is evident from (3) that o<u?<1; where necessary we shall take the positive 
root of u?, so that o<u< 1 also. 
* It will be noted that many of the arguments used here do not necessarily apply only to comets 


possessing particle swarm structure. However, this particular approximation cannot easily be 
made in other theories. 
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Substituting for r in (2) from (3) and putting s=tan 46, we obtain 


ht ylagt= [as (4) 


a (1 —u®)? 
where a=tan $a. 


In the case « = 1, for which u=o always, this integrates immediately to 
o(t—t')=2(s—a) + §(s*—a*) (5) 


where o =h/q? (independent of x). However, since we have chosen the origin of 
the time scale to be the instant of perihelion passage, t=o whens=o. Therefore, 
from (5), 


ot’ =2a + $a’. (6) 
Consequently, in the case « = 1, 
ot = 25 + ¥s%. (7) 


3. Locus at one instant of particles with different values of x.—Substituting for r 
in (3) from (1) for a particular value of «, we obtain 


u 1—«\'2 
3-Gza) (8) 
By substituting for ¢’ from (6) and for s from (8) in (4) and integrating, we arrive at 
the result 


ot=(1 +a*)(s—a) F + log a 
\ 


+ 
1—u® 2u I-u 

am 3 
, &—a) {+,- Zo i 

u2 I—u 2u I—u 

—ay? 

4 2a(s a) 

I—-u 
+ 2a+ $a’. (9) 


Since o<u<1 in (g) we can expand this expression in powers of u. It can be 
verified that if u is allowed to tend to zero the result (5) is obtained, corresponding 
to k=I1. 

4. Predicted appearance of a tail_—We could if we wished obtain explicitly an 
equation in polar coordinates of the locus of dust particles of varying « by sub- 

' stituting in (g) the expressions that we have for u and s in terms of r and 6 but 
independent of x. However, it has been found more convenient to continue using 
our present variables. 

To obtain a better indication of the distribution of dust particles after collisions 
within the head, it is necessary to consider what happens when such collisions 
take place at more than one instant. To do this it is only necessary to permit @ 
to vary in (9). 

Equation (g) is plotted in Fig. 1 in a representative case, using various values of 
t and assuming that the tail began to form at perihelion. For each value of t 
three values of a were chosen, corresponding to « equal to 0°, 15°, and 30°. Aset 
of values of u were then taken and (g) solved for sineach case. A short programme 


was written for the IBM 650 computer for this purpose, equation (9) being solved 
for s by iteration. 
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5. Some properties of a one-dimensional tail.—The calculation of the angle 
between a tail produced at a single collision point and the path of the cometary head 
can readily be made. From (3), we have 


x=rcos@=q(1—s*)/(1—u*) 
and 
y=rsin 0=2qs/(1—u?) 


whence, differentiating these with respect to s and dividing, we find 


dy _ 2sk+1—u? 
dx (1—s*)A—s(1—uw?) 














PRODUCTION OF TAIL 
| MATERIAL WITHIN THE HEAD 
CEASES AT THIS POINT 


PARABOLIC PATH OF COMET _ 


PERIHELION OF COME TARY 
ORBIT—FORMATION OF 


“eee eee = He 
q aoe TAIL BEGINS. 








Fic. 1.—The calculated distribution of dust particles at varicus times ot after perihelion passage. 
It is supposed that the tail begins to form at perthelion, and the dotted lines show later distributions of 
dust particles emitted from the head of the comet at this instant. The dashed lines shew the distri- 
bution of the particles emittea at perihelion plus 30°, while the thick lines indicate the distribution of 
particles emitted at perihelion plus 15°. If collisions are considered to be taking place almost contin- 
uously between the extreme positions, most of the tail material can be expected to be concentrated 
about the 15° loci at any later time. Each of the lines represents a graph of equation (9) for fixed 
t and a and can be regarded as extending to infinity, although that part of each line further from the 
Sun than about 2°5q is not shown. The arrows show the direction in each locus of inc reasing 
u(—>1) and decreasing n(—— ©). 
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where A=u du/ds (derived from (g)). Putting u=o, we have that the inclination 
of the tail at the head of the comet is 
25Aq + I 
(1—s*)Ay—s (11) 
where A, is the value of A at u=o. The inclination of the path at the point 
(r,@) is—1/s, so the tangent of the angle required is Ap/(s\g +1). To find Ay we 
differentiate (g) with respect tos and putu=o. This reduces to 
15 (1+5°) 
2(s—a){10 + a? + 3as + 6s*} 
and the expression for the tangent of the angle between the tail and the direction of 
motion of the head is obtained in the form 
15 (1 +s?) (12) 

353 + 6as? + (4a? — 5)s +20a+ 2a" 

In any particular case this will give the angle at which the tail leaves the head, 
measured from the direction of motion of the head, at any time after the comet 
passes a single collision point for values of s>a. The required value of s corres- 
ponding to a particular time ¢ (measured, of course, from perihelion passage) is 
found by solving (7) fors. As an example we can take the simple case a=o, that 
with the collision point at perihelion. Since the angle between the tail and the 
direction of motion of the head is a right angle when the cubic denominator 
vanishes, we immediately see that in this case this is so for s=(5/3)"*, which 
corresponds to an angle 6 of roughly 1044°. Equation (7) gives the value of ot as 
(28/9)(5/3)!? or about 4:0; up to this time the cubic denominator is negative, so 
the angleis obtuse. For all ot greater than this value the angle is acute, and tends 
to zero as s-> 0, that is as the comet moves away from the Sun towards aphelion. 
This is of course an expected result and accords with observation. 

To find the behaviour of the locus as r->oo we merely have to investigate 
equation (9) for u>1. For simplicity we again take the special case a=o, but 
similar arguments are applicable in the general case. Writing 


tang= 





y= - 





I I I+t 
so— 14 eee 
I—u 2u I—-u 


and 

I I I I+u 

B= —(——,, — — log —_ 

er au Tu 
ot=sA+s°B. (13) 
It can be easily verified that A > B > 0 for all uin the permissible range, and that 
as u tends to unity both A and B diverge. However, since s, A, and B are all 
positive and at any instant ot is constant for all u, sA must be less than ot for all u. 
Thus s must tend to zero in order to preserve this inequality. Consequently, since 


A>B>o, s°B tends to zero. By a similar argument it may be shown that the 
logarithmic part of sA tends to zero also, and so 


equation (g) reduces to 


s/\(1—u?)>ot as ut. 
Now (3) may be rewritten in the form 
r/(1+s?)=q/(1—u?) 
from which follows 
y=rsin@ =29s/(1 —u?). 
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‘This measures the distance y of any point on the locus from the initial line. Thus 
asr—> x 


y—>2got=2Vt 


where V is the velocity of the comet itself at perihelion. However, since the 
particle for which « =0 (i.e. that for which the gravitational force and the radiation 
pressure exactly cancel) travels in a straight line perpendicular to the initial line 
with constant velocity V, we see that the end of such a tail moves away from the 
axis with twice the velocity of the free particle. 

From (10) it is easily demonstrated that the inclination of the tail is positive for 
«<0. The conclusion is that the tail at infinity is parallel to the initial line. 

Similarly, in the general case the end of the tail moves away from the line 
joining the Sun and the collision point twice as fast as does the free particle. In 
fact, asr-> x 

y' =rsin(@—«)>2V'(t-t’)(1+ a?) 


where V’ is the velocity of the comet at the collision point. V'(1 +a?)~"? is of 
course the component of V’ perpendicular to the radius vector. 

6. Applications to Comet Arend—Roland.—Some of the preceding results have 
particular application to Comet Arend—Roland 1956h. It has been pointed out 
(2, 3) that at the time of the appearance of the anomalous sunward tail the Earth 
was passing through the plane of the cometary orbit. According to the elements 
compiled by Candy (4), the Earth was at the ascending node of the orbit on 1957 
April 25, at an angular distance of about 51°15’ past the cometary perihelion. 
Taking the Earth’s distance from the Sun at this time to be one astronomical unit, 
a rough calculation shows that the comet itself was then at an angular distance of 
about 82 past perihelion. 

By ditferentiating (11) with respect to a, and making use of the expression we 
have for Ay, it is easy to show that at any particular point of the orbit defined by 
keeping s constant the inclination tan ¢ of a one-dimensional tail at the head of a 
comet is a monotonically increasing function of a. Now the gradient of the line 
joining the Earth and the comet on 1957 April 25 was about 0-42. Putting 
s=tan41 and tan ¢6=0-42 in (11) gives a value of a of about 0-27, which corres- 
ponds to a negative value of « of about —30. Since tan ¢ is a monotonically 
increasing function of a, this means that if the tail began to form at an earlier stage 
than 30 before perihelion passage, the particles emitted from the head at the 
earliest collisions points will lead to values of tan ¢ smaller than 0-42. The 
situation is illustrated in Fig. 2, which demonstrates clearly that, at the time that 
the Earth passed through the plane of the orbit, that part of the tail which had been 
produced by collisions occurring before the comet reached an angular distance of 
30° before perihelion would have appeared from the Earth to be on the sunward side 
of the comet’s head. ‘That part produced by later collisions would appear as a 
normal tail. ‘The similarity between this diagram and Armstrong’s conjectural 
figure (3) is clear. 

The work of Lyttleton has shown that it is not unreasonable to expect some 
collisions to take place so early (1, p. 117), although most collisions can be 
expected close to perihelion, which explains the simultaneous appearance of a 
normal tail. It can also be seen that, since the Earth was at the ascending node 
of the cometary orbit on 1957 April 25, the ‘spike’ would have appeared to point 
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slightly to the south of the Sun just before the node was reached, and slightly to 

the north just afterwards. This is again in accordance with observation (2). 
At first sight it might be expected that, as in the case of the rings of Saturn, the 

anomalous tail would be too narrow to be visible to an observer in its own plane. 





Xe 


EARTH 


LOCUS OF PARTICLES 
f— EMITTED AT 


« = ~40° 


t 


LOCUS OF PARTICLES 
EMITTED AT 


« = -60° 


COME TARY OE ta.u 
Cant <> ©.3167 a.u. 


51°15! 
= 62° 


A PERIHELION 











Fic. 2.—The situation as the Earth passed through the plane of the orbit of Comet Arend- 
Roland on 1957 April 25. The tail as given by the IBM 650 programme (see Section 4) is shown 
as in Fig 1, together with two loci of particles emitted before the comet reached 30° before perihelion. 
The figure demonstrates how such particles appear on the sunward side of the head, while particles 
emitted later go to form the normal tail. The figure by Armstrong referred to in the text (Section 6) 
ts reproduced below for comparison. 





Plane of Ecliptic \ 


However, we have throughout taken the probiem to be two-dimensional to bring 
it within the range of analysis. Considerable dispersion of the tail outside the 
orbital plane (but with the greatest density at the plane) is to be expected, since the 
cometary head is not the point-particle we have assumed it to be, but a swarm of 
particles, with many of these necessarily having small motions perpendicular to 
the plane of general motion. Also, we have not taken into consideration any 
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additional components of velocity of the ejected particles perpendicular to the 
orbital plane, whereas with random collisions we must expect such components to 
add to the dispersion of the tail. 

Nevertheless, a cometary head can subtend only a small angle at the sun, so 
the consequent dispersion perpendicular to the plane is not likely to be large 
compared with the dispersion in the orbital plane revealed by our analysis. In 
addition, since only a small proportion of the tail can be expected to be formed 
earlier than 30 before perihelion, it is not unlikely that, to an observer in the 
orbital plane, only those particles in or near the plane will be dense enough to give 
a visible effect. ‘hus the narrowness of the spike in contrast to the dispersion of 
the tail should cause little surprise. 

The author is greatly indebted to Professor H. Bondi and to Dr R. A. Lyttleton 
for much advice and encouragement in the preparation of this paper, and to Mr. J. 
K. Iliffe and the Applied Science Department of IBM United Kingdom Limited 
for generously providing computing facilities. ‘The figure from Dr Armstrong’s 
paper (3) appears by courtesy of Dr Armstrong and the Editors of Observatory, 
to whom also the author wishes to express his thanks. 


Since this paper was prepared, a paper by Osterbrock (5§) on the tails of the 
comets Baade (1954 h) and Haro-Chavira (1954k) has appeared. Calculations are 
shown which demonstrate that the tails of these comets, although lying in the 
orbital plane, did not lie in the radial direction from the Sun, and explanations are 
suggested in terms of the composition of the tail material. It seems that the 
present paper has some bearing on this problem, since the preceding analysis 
shows clearly that a cometary tail does not necessarily lie in the radial direction. 
King’s College, London: 


1958 October 22. 
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COLLISION BROADENING AND SHIFT IN THE RESONANCE LINE 
OF CALCIUM 


W. R. Hindmarsh 


(Communicated by the Director, Wniversity Observatory, Oxford) 
(Received #958 November 21) 


Summary 


The collision shift and broadening of the line Ca 4227 A due to an external 
pressure of helium have been measured. The line was formed in absorption 
by passing white light through calcium vapour in the presence of helium at 
various pressures of less than 1 atmosphere. The spectral resolution was pro- 
vided by a Babcock grating used in the new Oxford solar spectrograph. The 
half-intensity damping width of the line was found to be 1-72 + 0°12 x 107*°cm=! 
per atom per cc of helium, and the shift 0-05 +0°04 x 10-*° cm™ per 
atom per cc towards the violet. The ratio of broadening to shift is therefore 
>19, compared with the value of 2°76 predicted by the Lindholm theory. 
The shift is also predicted to be towards the red. Possible explanations of the 
discrepancy are discussed and it is shown that the short-range repulsive 
forces between calcium and helium atoms must be involved as well as the long- 
range van der Waals forces. With the use of Lindholm’s theory, coefficients 
for both types of force are evaluated from the measurements. Finally it is 
shown that the solar damping coefficient for Ca 4227 A is <2°2 x 107° cm™ 
at t=0'1, and that the collision shift makes a negligible contribution to the 
solar red shift. 





Introduction.—Collision broadening is an important process in the formation 
of stellar absorption lines, and a knowledge of the damping coefficients of such 
lines would remove an important unknown from the theory of line contours. 
Collision shifts have been discussed as a possible explanation of the fraction of 
the red shift of solar lines which is unaccounted for by the Einstein gravitational 
shift (1, 13). ‘The work described in this and the following paper forms the 
first two of a series of investigations into collision effects in atomic spectra. 
Stellar absorption lines are formed, in general, in atmospheres of rather low pressure, 
and the present investigations are therefore aimed at the pressure region below 
1 atmosphere. In stellar atmospheres, collision broadening is probably due in 
many cases to collisions of the radiating atoms with both neutral and charged 
particles. ‘The experiments here described were designed to measure neutral 
particle broadening and shift only. 

The most useful theory of collision effects in spectra developed to date is 
that due to Lindholm (9g, 10). A similar theory was given independently by 
Foley (4). This theory gives the following expression for the absorption 
coefficient «(v) of an atom absorbing in the presence of N perturbing atoms 
per cc:— 
tee constant 

= BP 
8 and yare constants, and visincm~!. This expression represents a symmetrical 
contour of damping shape and half intensity width 2y shifted with respect to the 
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line centre at v= v, by an amount 8. Assuming a van der Waals interaction—an 
attractive force whose potential is proportional to the inverse sixth power of 
distance between the radiating and perturbing atoms—the values of 8 and y are 
B= — 4:90 x 10-7C,25035N 
y = 6°75 x 10-7C,7 55 5N ; 
v is the mean relative velocity of the colliding and radiating atoms; C, is the 
interaction constant defined by v= — C,/r®, dv being the change in wave number 
of the radiation caused by a perturbing atom at distance 7. Lindholm (9, 10) 
and also Unséld (15) have discussed the range of validity of this expression. 
They find that the approximations involved are good for those parts of the 
line contour for which 
pos 
(v— v9) < 6-09 x 10714 rans 
6 
Substituting for C,, this inequality becomes 
‘ N\ 12 
(v — v9) < 5:00 x ro-re92 (™) R (1) 
Y 
y/N frequently lies in the range 10-!*—10-*%°cm~! per atom per cc, and 
v~2»x10°cemsec™! for most laboratory experiments with light gases. In these 
circumstances the condition is (v—vy)<45cm~, or about 8 A at 4227 A. The 
low pressure formula for «(v) may therefore only be applied to the comparatively 
narrow lines formed at particle densities of less than about 2 x 10!* per cc, for 
which y < 45 cm™!. 

Little previous experimental work on collision broadening and shift has been 
carried out at these low pressures. Kleman and Lindholm (6) have measured 
the broadening and shift of the sodium D lines produced by argon at low 
pressures, and found that the ratio 2y/8 was about — 4-0 as compared with the 
theoretical value of — 2-76. Kusch (7) has recently measured some iron lines 
broadened by atomic hydrogen in conditions which just satisfy the criterion, and 
finds that the ratio 2y/8 varies from —2-5 to — 3:6 for different lines. Kusch, 
however, claims an accuracy of only + 30 per cent for the ratio 2y/8 and + 40 per 
cent for his values of y/N. He has also measured the broadening only of the 
resonance line of calcium, A 4227, in an atmosphere of molecular hydrogen. 
Lang (8) has measured the broadening only of some neon and argon lines excited 
in constricted glow discharges at pressures below 2 cm mercury. 

It seems, therefore, that accurate measurements of collision broadening and 
shift would be of considerable value, both for the adequate testing of the 
Lindholm-Foley theory and for the elucidation of certain astrophysical problems. 
There are at least two possible methods of making such measurements. First, 
the line may be measured in absorption (with the use of a simple absorption tube 
for resonance lines of volatile elements, or a King furnace for other lines). This 
has the advantages that the damping coefficient of an absorption line is com- 
paratively easy to determine and that the line is produced in well defined 
conditions, and in the virtual absence of charged particles. The disadvantages 
are that, even with the best available diffraction gratings, pressures of a 
substantial fraction of 1 atmosphere are necessary to obtain accurate measurements, 
and that lines of excitation potential greater than about 10 volts cannot be 
obtained in a King furnace. The second possible method is to excite the line in 
emission in a gas discharge and use interferometric techniques for the 
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measurements. With this method, very low pressures may be used, though in that 
event the broadening is exceedingly difficult to measure. Further, the presence of 
charged particles in the discharge makes necessary some investigation of their 
influence on the line. This paper provides an example of the former type of 
experiment, and the following paper an example of the latter. 

Experimental method.—The resonance line of calcium, A 4227 A, formed in 
absorption in an atmosphere of helium, was examined, and the collision shift 
and broadening determined. The calcium and helium were contained in a 
steel absorption tube 15cm long and 2cm in diameter, which was heated to a 
temperature of about 400 °C by means of an electric furnace. The temperature 
for a given heating current was measured in a separate experiment when a 
(platinum) — (platinum-rhodium) thermocouple was inserted in the tube. This 
could not be done in the presence of calcium vapour, as the latter reacted 
chemically with the thermocouple wires. At a temperature of 400 °C the vapour 
pressure of the calcium was sufficient to give an absorption line of fairly small 
but easily measurable central intensity. The ends of the tube were cooled by 
means of water jackets, and apiezon wax was used to attach glass windows. The 
influence of the non-uniform temperature in the absorption tube on the line 
contours was estimated. Errors in the measured damping coefficient due to this 
cause are not more than 3 percent. Physically this small effect arises because of 
the very rapid variation of vapour pressure with temperature, so that only a thin 
layer in the hottest part of the tube is effective in forming the absorption line. 
The tube was connected to a mercury manometer and to a cylinder containing 
helium through a liquid-oxygen cooled charcoal trap. The whole system could 
be evacuated by means of a rotary-oil pump. The purpose of the charcoal trap 
was to ensure that no trace of oxygen entered the absorption tube, as calcium 
reacts very rapidly with it to form calcium oxide, a stable and refractory substance. 

A tungsten strip-filament lamp provided a source of continuous radiation. 
The light passed through the absorption tube in a collimated beam, then through 
a combination of colour filters to isolate the spectral region required, and an 
image of the strip was formed in a plane in which a rhodium step filter could be 
placed for calibration purposes. Finally, the filament was re-imaged on the slit 
of a plane-grating Littrow spectrograph of 12 m focal length, the new Oxford 
solar spectrograph. ‘The grating was ruled by H. W. Babcock. It has 600 lines 
per mm and an aperture of 20cm. It was used in the sixth order, where the 
blaze was effective, and the dispersion of about 7mm per A was measured from 
iron-are spectra taken for the purpose. The spectrograph slit was covered with 
a diaphragm which divided the slit length into four sections, each of which could 
be opened or closed by means of a shutter. On the lower two of these divisions 
fell the image of the strip filament lamp, while on the upper pair light from a 
G.E.C. cadmium lamp could be directed by means of a right angled prism. The 
cadmium green line A 5086 A in the fifth order fell about 5 cm from the calcium 
line A 4227 A in the sixth order, and the long-wavelength hyperfine-structure 
component of this cadmium line was used as a reference line to evaluate spurious 
shifts due to changing atmospheric pressure and other causes. 

The experimental procedure was as follows. Helium was admitted to the 
tube at one of the three pressures 6cm, 40cm, or 75cm of mercury, the tube 
already containing calcium. The tube was maintained at a temperature of about 
400 °C. The continuous spectrum of a strip-filament lamp was then exposed 
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simultaneously with the cadmium green line, each through one of the four 
divisions of the slit length. ‘Then the shutters were moved to open the other 
two divisions, the helium pressure changed to one of the other two values already 
mentioned, and a similar exposure taken. One of the two helium pressures was 
always 6cm of mercury. Exposures for plate calibration were made in the 
usual way, with the use of a calibrated rhodium step filter (3). With Ilford 
Zenith plates a time of 2 hours was required for each exposure when the slit 
width was o-12mm. 5 plates were obtained for the measurement of the shift 
between 6cm and 40 cm mercury pressure, and 5 for the shift between 6 cm and 
75cm. 6, 5 and 4 spectra were used to evaluate the broadening at the pressures 
6 cm, 40 cm and 75 cm respectively. 

Measurement of the shift. ‘To evaluate the wavelength shift of the calcium 
line caused by changing the pressure of helium in the absorption tube, the 
positions of the calcium and cadmium lines were measured on a modified Hilger 
measuring micrometer whose plateholder could be moved in two mutually 
perpendicular directions. ‘The appearance of a plate is shown diagrammatically 
in Fig. 1, where the shifts have been exaggerated for the sake of clarity. The 














Ca ltirst exposure) | 
Ca (second exposure) | ———+ | 
| —|caltirst exposure) 
| ee ca(second exposure) 
Fic. 1.—Diagrammatic illustration of the appearance of a plate. 


shift in mm between the two cadmium exposures was subtracted from the shift 
in mm between the two calcium exposures to give the true collision shift of the 
calcium line, which was converted into cm~! and A with the use of the measured 
value of the dispersion. Mechanical displacements of the plate, slit or grating 
during an exposure are clearly eliminated by this procedure. Consider now the 
effect of changing atmospheric density. If the latter changes by Ap then the 
changes in wavelength of the calcium line and of the cadmium line are given by 
* = =(",.—1) 1 and Aci =(p.4-1) ap ; 
A a Pp ay d 
where yz... and j,., represent the refractive indices of air at A... and A,., respectively. 


Then 





Ca Mca! 


A 
Ad,,. = AA,, — .§ —— 
- Me Aca Mca! 


and if the dispersions of the spectrograph at the two wavelengths are D,,, and 
D,., mm/A and the shifts in mm on the plate are Ax,., and Ax, 
x Nea (Hea — I) Des 
- d eS 

me a (Mca — 1) Dew 
= Ax,,, to an accuracy of 1 per cent. 


Ax,.. = 
a 
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Hence the spurious shifts due to changes in air density are eliminated by subtract- 
ing the cadmium shift in mm from that of the calcium shift in mm. Further 
errors may arise because of lack of perpendicularity of the two stage motions 
and inaccuracy in the alignment of the plate on the microscope. These were 
eliminated by ensuring that the distance through which the microscope was 
moved parallel to the length of the spectral lines (marked by the square brackets 
in Fig. 1) was the same for both the cadmium and the calcium shifts. An estimate 
was made of the influence of curvature of the spectral lines, and this was found 
to be negligible. Each plate was measured twice, and inverted between the two 
measures. ‘There was no significant difference between these two measures. 

Table I shows the results obtained. The quoted error is in each case the 
standard deviation of the mean. The pressure has been converted into particle 
density with the use of the measured temperature, which was 400 °C for the 
6 cm pressure exposures, 392 °C for the 40 cm and 416 °C for the 75 cm. 


TABLE I 


Shift between 0-86 x 10!8/cc Shift between 0°86 x 1018/cc 
and 5-80 x 10!8/ce atoms of helium and 10-48 x 1018/cc atoms of helium 
1607* an“ mA 10-*cm7! mA 
+15°1 —2°7 + 10°5 —1'9 
+ 5°6 —1°0 + O°5 —o'l 
—24°0 +4°3 + 73°: 13 
+ 73 —13 +15°1 —2°'0 
= ea +04 + 3°9 ~O'7 


Mean+ 0°4+3°4 —o1+06 Mean+ 7°4+2°3 —1'2+04 


From these results we find a mean value of the quantity B/N, assuming that 
the shift between x atoms of helium per cc and y atoms of helium per cc is equal 


to that between x—y atoms per cc and zero atoms per cc, provided y<x. The 
result is 


B/N =0-05 + 0°04 x 10-*° cm~! cm, 

Measurement of the damping coefficient—The calcium absorption lines 
produced at the three pressures of helium were measured with the use of a 
Moll-type microphotometer (12). The characteristic curve of each plate was 
plotted from microphotometer tracings of the step-wedge exposures. The curves 
for all plates but one were averaged together. The curve of the remaining plate 
was used for that plate, since it was significantly different from the others. Two 
tracings of each calcium exposure were taken from portions of the absorption 
line, and photographic density reduced to intensity of light in the usual way. 

It was assumed that the line contours could be expressed in the form 
I(v)=1,e-*” where «(v) is the atomic absorption coefficient of calcium modified 
by the apparatus function of the spectrograph and multiplied by the number of 
atoms per cm? of the absorption tube cross section. a(v) is the optical depth 
of the calcium vapour at frequency v, and it was computed from the observed 
I(v). The resulting curve was fitted to a Voigt profile at the o-5 and o-1 intensity 
widths with the use of the tables given by van de Hulst and Reesinck (§), and 
the half intensity widths of the damping and Gaussian components determined. 
The fit to a Voigt profile was in all cases better than 3 per cent of the maximum 
value of «(v). A typical curve for each of the pressures is given in Fig. 2, and 
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Fic. 2.—Observed values of a(v), together with the fitted Voigt profiles, for a typical plate 
for each pressure. (a) 6 cm helium pressure; (b) 40 cm helium pressure; (c) 75 cm helium pressure. 
Ordinates are on an arbitrary scale, and abscissae are in cm™' measured from the line centre. 


points from both the tracings for the plate chosen—that at the head of its respective 
column in Table Il—plotted. Table II shows the individual plate measures, 
and the errors are standard deviations of the mean. There are no significant 
differences among the Gaussian half-intensity widths measured at the three 
different pressures. <A plot of the damping half-intensity width against particle 
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Tasie II 


6 cm helium pressure 40 cm helium pressure 75 cm helium pressure 
0°86 x 10'* atoms helium/cc 5°80 x 1018 atoms helium/cc 10-48 x 10'* atoms helium/cc 


Damping Gaussian Damping Gaussian Damping Gaussian 

width width width width width width 
on cm-! cm-! cm} cm! cm~? 
0060 0°148 O°15§2 0°147 0°216 0°202 
0089 0'142 O°175 O°195 0°225 0°173 
0076 o'102 0°156 0°093 o'241 O°175 
0083 0°129 0°138 0099 0'271 0°065 
0048 0°128 0'138 O°1l4 

0086 0°094 


©°073+0°006 0°124+0°008 0°152+0°006 0°130+0°017 0°238+0°0I10 0°154+0°026 





03 T 





0.0 | } i i 
0.0 2.0 4.0 6.0 8.0 . 12.0 








Fic. 3.—Plot of damping half-intensity width against helium particle density. The straight 
line is the least squares line through the three observed points. Ordinates are in cm and abscissae 
in units of 10'* atoms helium per cc. 


density of helium is shown in Fig. 3. The equation of the least squares line 
through the three points is 

Av=1'72 x 10-°N +.0°056. 
There is every reason to believe that the collision damping coefficient is zero at 
zero particle-density of helium, for ‘‘ self ’’ broadening of the calcium must be 
negligible at the very small vapour pressures of calcium involved in this 
experiment. The residual damping width of o-056cm~! at N=o is therefore 
attributed to the apparatus function of the spectrograph, and may be subtracted 
from the observed damping width to give the true collision width, so that the 
result is 

2y/N= 1-72 + 0°12 x 107" cm—! cm, 
The mean value of all measures of the Gaussian half-intensity width was 
Av, =0°134 + 0°010cm~—1, 
This is due to the Doppler width of the absorption line and the apparatus function 
2 
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of the spectrograph. As a check on the plausibility of this figure and the figure 
used above as the damping component of the apparatus function, the data given 
by Pierce (11) for the apparatus function of a Babcock grating were assumed for 
the grating used in the present experiment, and the Gaussian and damping 
components to be expected here were computed. ‘The expected damping half 
width was 0-046cm~!, compared with that observed of o-056cm~', while the 
Gaussian apparatus half width added to the calculated Doppler width was found 
to be o:166cm~', compared with the observed value of 0o-134cm~!. The 
agreement between the expected and observed values is as close as can be hoped 
for with the assumptions involved. 

Discussion of the results.—The simplest test of the Lindholm theory is to 
compare the experimental and theoretical values of 2y/8, the ratio of broadening 
to shift. The value obtained from the above results is 


19<2y/B< « 
while the theoretical value is 2-76. Furthermore, the observed shift is in the 
opposite direction to that predicted. There is therefore a marked disagreement 
between theory and experiment. 

As a further comparison of theory and experiment, an attempt has been 
made to evaluate the van der Waals attraction between calcium and helium 
atoms, and so to estimate a theoretical value of 2y. The expression for C, due to 
London was used in the form given by Unsold (equations 82, 49 of 15). The 
results of Hartree-Fock calculations of the oscillator-strengths were used where 
these were available (Trefftz, 14). The remaining oscillator-strengths were 
calculated by the method of Bates and Damgaard (2). The result is 

2y/N=1-74 x 10-* cm-! cm? 
so that the agreement here is good. 

Before seeking explanations of the discrepancy between theory and 
experiment for the ratio 2y/8, it is necessary to investigate whether the low 
pressure approximation ts valid in the conditions of the experiment. Substitution 
in equation (1) gives the condition for validity as 


v—vy<48cm~, 


Since the values of «(v) in these experiments were in all cases indetectably small 


1 the criterion is well satisfied. 

There are two possible explanations for the major discrepancy revealed by 
these measurements. First, the Lindholm theory may be incorrect even at low 
pressures. ‘lhe theory treats the collision process classically. Lindholm has, 
however, considered the modifications introduced by a quantum mechanical 
treatment of the collisions in the special case of the sodium D lines broadened 
by argon, and has found that only quite small corrections to the classical 
expressions are needed. It seems unlikely that the case of the calcium resonance 
lines broadened by helium will be very different, although this explanation of the 
discrepancy cannot be wholly discounted. The second possible cause of the 
disagreement is that the interaction of the calcium and helium atoms may not be 
described by a van der Waals (r~*) law. ‘To investigate this, it is convenient to 
evaluate from the experimental value of 2y/N the collision cross-section o, for 
this gives the order of magnitude of the distance of closest approach of the 


for (v—v,)>1cm 
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calcium and helium atoms. Now y/N=vo/c, from which we find 

ao = 1°3 x 1075 cm’, 
and the collision radius is 2-1x 10-*cm. Since the radius of the helium atom 
is about 1-o A and that of the calcium atom about 2:3 A, the collisions effective 
in broadening the line are rather close collisions, and at these small interatomic 
distances the repulsive forces between the atoms must be operative. Thus the 
assumption of an r~® interaction is probably invalid. The repulsive forces 
contribute to the broadening, and produce a shift in the opposite direction to 
the attractive r-® forces. It is clear, therefore, that the observations can be 
accounted for qualitatively in this way. More quantitatively, we may carry 
through the Lindholm theory for an interaction of the form 
Co, Gis 
fF = 
An r~™ form for the repulsive forces is frequently adopted in molecular theory. 
When the Lindholm theory is developed for this interaction, we find 


év=— 


_ "75 x 10>* Cs /5g3/5 + 1°94 x 107% Cic* 1199 1l 


Y 
N 
B 


NW 499% 1077 C25 v55 +. 0°57 x 10-9 Cy Fy 11, 


Substitution of the experimental values for y and 8 then yields 

Cy =5°31 x 10-* cm-! cm? 

Cig = 1°20 x 10-* cm-"cm® 
compared with the calculated value of C, (corresponding to the calculated 2y/N 
given above) 

Cy =2°13 x 10“ cm! cm*, 
There appears to be no simple method of calculating C,,. There is not now 
good agreement between the calculated and measured values of C,. The position 
is that we may retain the Lindholm theory by modifying the force law between 
calcium and helium atoms and accepting the disagreement with the calculated 
value of Cy, or we may retain a pure r~® law (though the measured collision 
cross-section makes this procedure of very dubious validity) with the good 
agreement between calculated and observed values of Cy, and suppose that the 
Lindholm theory is correct for the broadening but incorrect for the shift of the 
line. In either case the present situation is unsatisfactory, and can probably 
only be improved by making an extensive series of measurements on a variety 
of spectral lines. 

It remains now to discuss the astrophysical significance of the results. In 
stellar atmospheres, the collision effects in Ca 4227 A are due almost entirely to 
hydrogen atoms and electrons, while the measurements described in this paper 
refer to helium atoms. If we assume that the force between the atoms is of the 
van der Waals type (r~*) for both helium and hydrogen, then the ratio C, 
(hydrogen)/C, (helium) has been calculated by Unsdld (15) to be approximately 
30. Taking a layer in the solar atmosphere where the density of hydrogen is 
5 x 10'6 atoms per cc (optical depth s=0-1) and the temperature 5000 °K, then 
using the measured value of y (helium), we find 


y(hydrogen) = 1-9 x 10-*cm-1" 
2* 





20 Collision broadening and shift Vol. 119 


This computation takes no account of the probable importance, revealed by the 
measurements, of the repulsive forces between calcium and helium atoms. 
There is no simple method of comparing these forces for hydrogen and helium. 
However, since the range of the repulsive forces in the hydrogen case must be 
smaller (because of the smaller radius of the hydrogen atom) while the strength 
of the attractive van der Waals forces must be greater (because of the larger 
polarizability of hydrogen) the repulsive forces are expected to be relatively 
less important in the hydrogen case. Thus the value for y(hydrogen) is an upper 
limit. The contribution of the electrons may be estimated from the known Stark 
coefficient of Ca4227 A, taking the electron density in the layer concerned 
to be 6x10!"%cm-* and using the Lindholm theory. The result is 
y(electron)=2-8x 10°-4cm~!. An upper limit for the total damping coefficient 
in this particular layer of the solar atmosphere is then given by 
y(Sun) = 2-2 x 10-3em=}, 
An upper limit for the solar collision shift may be found by using the measured 
y(helium), calculating S(helium) from the Lindholm theory, converting to the 
conditions of the solar atmosphere as described above, and adding a calculated 
contribution for the electrons. The result is 
B(Sun)=1-7 x 10-3cm-!. 

If the measured shift is taken for 8(helium), then a much lower value for 8(Sun) 
isobtained. Inany event, the collision shift for Ca 4227 Aisa negligible component 
of the solar red shift, which is of order 5 x 10-*cm~', and cannot account for the 
difference between observed and predicted solar red shifts. 
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Summary 


The collision shifts of two argon lines, A 7067 A and A 6965 A, have been 
measured, the lines being excited in a D.C. glow discharge in pure argon. 
The lines were examined by means of a Fabry-Perot interferometer crossed 
with a Hilger prism spectrograph. The interference ring diameters were 
measured for a variety of argon pressures and discharge currents. It is 
shown that the measured shifts were due to collisions between two argon 
atoms, and not to the influence of charged particles. The shift for A 7067 
is —0o-61+0°05 x 10-*° cm=! per atom per cc of argon, and for A 6965 is 
—o'65 +0°'05 xX 107° cm! per atom per cc of argon, where the negative 
sign indicates a red shift. "These results are shown to imply that the short 
range repulsive forces between argon atoms can be of little importance in 
producing the shifts, in contrast to the case of the calcium line A 4227 A 
perturbed by helium. 





Introduction.—One of the two possible methods of investigating collision 
effects in atomic spectra discussed in the preceding paper is to excite the spectrum 
in a gaseous discharge, and examine the spectral lines with the use of an inter- 
ferometet. Such a method is the only one available for a spectrum of high 
excitation potential, as is the argon spectrum, where thermal excitation is 
impossible. If the results obtained in such measurements are to be physically 
interpretable, it is necessary to establish whether collision effects are due to 
neutral particles, electrons or ions, or all three. This paper describes the 
measurement and interpretation of collision shifts for the two infra red argon 
lines A 7067 A and A 6965 A at pressures of argon below 3 cms mercury. The 
transitions giving rise to the lines are 

3P 45(45 1°)—3p 4P(4P 8), A 7067 
and 3P 45(45 1°)—3P 4P(4P 9), A 6965, 
where the notation is that of Revised Multiplet Tables (3). These lines should 
not suffer from “ resonance” broadening to any appreciable extent, but, since 
the lower level is metastable, they may suffer from self absorption broadening. 
Neither resonance nor self absorption broadening can produce any shift in the 
conditions of this experiment. 

Experimental method.—The discharge tube is shown in Fig. 1. It was made 
of ‘ Pyrex’ glass, and was 1 cm in diameter. The anode was a cylindrical 
copper insert (C) shaped to fit the ground glass sockets (G). Sealing wax 
was used to make the joints vacuum tight. The cathode consisted of a copper 
block (B) inserted into the copper-glass seal (S). The window (W), attached 
to the tube by means of black wax, allowed light from only the positive column 
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of the discharge to emerge. ‘The light from the cathode region, where the 
macroscopic electric fields are high, did not enter the optical system. The 
discharge tube was cooled by means of a water bath maintained at about 12 °C. 

The discharge tube was attached to a vacuum line, and the argon pressures 
measured by means of a manometer containing diffusion pump oil of density 
1:oggmcm *. A liquid oxygen cooled trap lay between the manometer and 
the discharge tube. 

The voltage developed across the discharge was about 300 volts, and was 
almost independent of current and pressure in the range 50-150 mA and 0:8-23°5 
mm of mercury. 

Light from the positive column was passed through the interferometer in a 
collimated beam, and the interference fringes were focused on the slit of a 
Hilger medium spectrograph with glass optics by means of a Zeiss achromat of 
focal length 84cm. Light from a G.E.C. cadmium lamp could also be passed 
through the interferometer, and the constancy of the interference-ring diameters 
of the cadmium red line A 6438 used as a check on the stability of the inter- 
ferometer. 

The interferometer plates were 6 cm in diameter and made of fused quartz. 
Silver films of estimated reflectivity 87 per cent and transmissivity 11 per cent 
at 7000 A were deposited on them by evaporation. These gave a theoretical 
instrumental half-intensity width of 0-005 cm™! with a 3 cm spacing between 
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the plates, and 0-004 cm™! with a 4cm spacing, compared with the Doppler 
width of the argon lines at 300°K of 0o-03cm~!. The room containing the 
interferometer and spectrograph was temperature controlled. The inter- 
ferometer was enclosed in an air-tight box to maintain the air density constant. 
Measurements were made with 3 and 4 cm spacers in the interferometer. 

Four exposures were made on each plate. First the cadmium spectrum was 
exposed. Secondly the argon spectrum was photographed at a certain discharge 
current and argon pressure. Thirdly the argon spectrum was again photo- 
graphed but at either a different current or a different pressure. Fourthly the 
cadmium spectrum was exposed again. Three argon pressures were used, 
namely o-8 mm, 12:0 mm and 23-5 mm of mercury, and three argon discharge 
currents, namely 50 mA, 100 mA, and 150 mA. Kodak I.R.E.R. plates were 
used, and exposure times for the argon lines were between 10 and 60 min, while 
for the cadmium red line the time was 3 min. 

Measurement of the shift.—On each plate the interference ring diameters 
of the two exposures of Cd 6438 were measured. If the interferometer 
displacement revealed by differences in the ring diameters between the two 
exposures amounted to more than 1/200 of an order of interference, the plate was 
rejected. On the plates which survived this test, the ring diameters of the two 
argon lines were measured for each argon exposure. ‘These measured diameters 
were corrected to allow for the very small displacement of the etalon during 
the argon exposures by assuming that the displacement expressed as a fraction 
of an order of interference deduced from the cadmium ring-diameters proceeded 
at a uniform rate. This correction was always much smaller than the measured 
argon shifts. The change in wavelength of the argon lines between two 
exposures on the same plate was then deduced from the corrected ring diameters 
by the method of rectangular array. When the two argon exposures concerned 
were on different plates, the shift between the corresponding cadmium exposures 
was subtracted from the corrected shift between the argon exposures to give 
the true shift of the argon lines. ‘This technique of measuring shifts is similar 
to that used in the measurement of isotope shifts in atomic spectra by Pery (2). 

‘Table I shows the results for changes of discharge currents at a fixed pressure, 
while Table II gives the shifts measured for changes of pressure at a fixed 
current. These tables show clearly that there is no measurable change of 
wavelength with discharge current, and that there are no significant differences 
between the results for different interferometer spacers. The error quoted 
is in all cases the standard deviation of the mean. 

Discussion of results.—The measured shifts may be due to changes in the 
static electric fields in the discharge, collisions between charged particles and 
argon atoms, or collisions between argon atoms. 

The estimated theoretical Stark effect of both lines is about 10-* cm~! for 
an electric field of 1 kV/cm, and dependence on the field strength is quadratic 
in form. It is most improbable that fields of this order can exist in the pvsitive 
column of a glow discharge, except in the space charge sheath formed near the 
walls. Since this sheath extends into the tube only for a distance of one mean 
free path, the light from it cannot make any observable contribution to the 
measured shifts. It is therefore concluded that Stark effects caused by the 
macroscopic electric fields of the discharge cannot be responsible for the observed 
shifts. 
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TaBLe I 


Change of wave number with current 


First and 
second 
currents mA 


Pressure 
mm Hg 


Wave- 
length 


150 
150 
150 
150 
150 
150 


spacer (no. of 
plates) 
10-* cm- 
+01 +0°4 (3) 
+02 +0°3 (3) 


+o'1 +01 (3) 
+0°2 +0°2 (3) 


+0°3+0°'1 (3) 
+0°3 +0°3 (3) 
+o'1 +0°3 (3) 

0:0 + 0°3 (3) 


+0°2 +0°2 (3) 
—0°3 +0°2 (3) 


+o°1 +0°3 (3) 
+0°2 +0°4 (3) 


TABLE II 


Shift with 3 cm Shift with 4 cm 


spacer (no. of 
plates) 
10~* cm-* 
—o'r+orr (3) 
00 +0°4 (3) 
—o'r torr (3) 
—0'3 +0°2 (3) 
—o'4 +0°4 (3) 
+o'2+0'1 (2) 


—o'r+0°5 (3) 

00 + 0°4 (3) 
—0°7 + 0°3 (3) 
—o°3 torr (3) 
—1ot+orr (3) 
+o'1 +0°4 (2) 
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Mean 


107° cm~! 
oo+o'2 
+o°r1 +0°2 
ooto'! 
o'o+0°2 
oo+to'2 
+o°2 +0°2 


070 +0°3 
ooto3 
—oa +o°2 
—-o3+0°!1 
—os os 
+0°'2+03 


Change of wave number with pressure (average of plates at 50, 100 and 150 mA) 


Wavelength 
pressures mm Hg 


plates) 
10~° cm! 


First and second Shift with 3 cm Shift with 4 cm 


spacer (no. of spacer (no. of 


plates) 
107° cm! 


—46+0°1 (6) 
—2°4+0°2 (6) 


—4'! +o'2 (5) 


—54+0°5 (5) 
—2°3+0°2 (6) 


—5'2+0°4 (4) 


Mean 


107? cm~? 
— 50 +0°3 
—2°4t0°2 


—47 +03 


—2°4+0°3 (6) —2'2+0°3 (6) —2°3t+0°2 


10 X 10-° cm~'=5-0 x 10-4 A at A=7000 A. 


Assuming the same value for the Stark coefficient of the argon lines, we find 
that Lindholm’s theory for broadening by electron impact predicts that an 
electron density of 10’ electrons per cc is required to produce a line shift of 
10° *cm™~}, if the electron temperature is taken as 3x 104°K. The value of 


the shift is very insensitive to the temperature assumed. Estimates based on 
measurements of the electron drift velocity in an argon discharge by Nielsen (1) 
suggest that the highest electron density likely to be encountered in our 
experiments is 5 x 10"! electrons per cc, so that the shifts due to electron impact 
are negligible. Effects due to ions will be very much smaller. (See, for 
example, Unsdld (4)). 

These theoretical estimates are supported by our observation that there is 
no detectable shift due to current changes. The conclusion is therefore drawn 
that the measured shifts are due to collisions between two argon atoms. The 
resonance forces cannot give rise to a shift, so that the van der Waals attractive 
and short-range repulsive forces must be responsible. As the theory predicts, 
the observed shifts are proportional to the pressure change at a fixed temperature . 
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The final result for the shift is 
A 7067 B/N = —0-61 +0°05 x 10°*® cm™ per atom per cc of argon 
\ 6965 B/N = —0°65 + 0°05 x 10°*° cm~! per atom per cc of argon. 

These values of B/N are substantially larger than that found for Ca 4227 per- 
turbed by helium, and, moreover, are in the expected direction, towards the red. 
Assuming that Lindholm’s theory is correct for the shift, and that the interatomic 
forces have a potential varying as the inverse sixth power of the distance between 
two argon atoms, we find that the collision radius necessary to give these shifts 
is 4x10 *cm. If the short-range repulsive forces are important, an even 
larger collision radius would be necessary to give the observed shift, since the 
repulsive forces tend to produce a violet shift. Since the radius of an argon 
atom in the ground state is only 1-56 10-8 cm, it seems probable that the 
short-range forces are not, in fact, very important in this case, in contrast to 
what was found for Ca 4227 in the preceding paper (5). 
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ANGULAR DIAMETER MEASUREMENTS OF THE RADIO SOURCES 
CYGNUS (1gN4A) AND CASSIOPEIA (23N5A) ON A WAVELENGTH 
OF 10°7 CM 
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Summary 


[he apparatus used in measuring the angular diameters of radio sources 
at a wavelength of 10°7 cm is described. Results for two directions at right 
angles are given for both Cassiopeia (23 N5A) and Cygnus (19N4A). These 
show that the former is circularly symmetrical and much the same size as it 
is at metre wavelengths, while the latter is markedly asymmetrical with a 
major axis that may be slightly larger at the higher frequency. ‘The position 
angle of this major axis has been found to be 109° +2°. 





Introduction.—The apparent angular diameters of the most intense radio 
sources have already been measured successfully at metre wavelengths (Smith 
1952, Mills 1952, 1953, Hanbury Brown et al. 1952, Jennison and Das Gupta 
1956). Subsequently the diameters of some of these sources have been measured 
at 500 Mc/s, (Conway 1956) and at 2800 Mc/s, (Rowson 1957). In both of 
these papers the high frequency results have been compared with the metre 
wave results and no very significant differences found. The present paper 
describes the apparatus used at 2800 Mc/s and gives the results obtained for 
Cygnus (1gN4A) (R.A. 19" 57™ 458, Dec. 40°35’) and Cassiopeia (23N5A) 
R.A. 23" 21™ 128, Dec. 58° 32’) for two directions at right angles. 

Apparatus.—The apparatus consisted of a two aerial interferometer working 
on a wavelength of 10-7 cm. Each aerial was a focal plane paraboloidal reflector 
illuminated by a dipole primary feed (Fry and Goward 1950) coupled to a wave- 
guide which ran along the axis of the pa:aboloid to a crystal mixer behind the 
reflecting surface. The larger reflector was 25 ft in diameter and consisted of 
half-inch wire mesh supported on a tubular iron frame giving a beam width of 
about 1° between half power points. The ro ft reflector was made of sheet 
iron and gave a beam width of 2}°. The aerials were supported on altitude- 
azimuth mountings, one being fixed to a concrete block while the other was 
mounted on a trailer, so that the length and direction of the baseline could be 
changed. 

A waveguide mixer and intermediate frequency preamplifier were mounted 
in a box behind each aerial, and the signal at intermediate frequency was fed 
along equal lengths of coaxial cable to the main amplifier. Here the signals were 
added together, further amplified, and the beat frequency between the two com- 
ponents obtained in a detector. 
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Separate klystron local oscillators were mounted on each aerial and about 
one per cent of the power from these oscillators was fed into the mixers via direc- 
tional couplers. The rest of the power, which would normally be dissipated in 
a dummy load, was fed into long lengths of coaxial cable to a crystal mixer placed 
between the two aerials, for baselines shorter than 800 wavelengths. For longer 
baselines it was found impossible to get sufficient signal through the cables 
because of their very high attenuation at a frequency of 2800 Mc/s. Thus 
for these longer baselines, the local oscillator at the distant aerial was radiated 
from a 4ft paraboloidal aerial and received on a similar aerial placed near the 
coaxial crystal mixer. This mixer produced the beat frequency between the 
two klystrons which was amplified by a selective amplifier at 20 kc/s and fed to 
a Foster—-Seeley discriminator, the d.c. output of which was used to control the 
frequency of one of the oscillators to keep the beat frequency stable to within 
+1 ke/s (Fig. 1). 
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Fic. 1.—Block diagram of apparatus. 


Interference fringes were obtained by selectively amplifying the intermediate 
frequency amplifier’s video output and detecting it with a coherent detector 
which used as its reference wave the 20 kc/s beat frequency between the two 
klystrons. The signal from this detector was amplified in a d.c. amplifier and 
displayed on a pen recorder. The total power from the intermediate frequency 
amplifier was also recorded for gain calibration. The theory of an interferometer 
using local oscillators of slightly different frequencies in the two channels has 
been given by Hanbury Brown, Palmer and Thompson (1955). 





28 B. Rowson, Angular diameter measurements of Vol. 119 


When using East-West baselines longer than 400 wavelengths it was found 
that the fringe frequency was inconveniently high, while for observations at 
transit with North-South baselines the fringe frequency was zero. To over- 
come these difficulties a system known as the ‘phase shifter’ was developed, 
which changed the frequency of the reference signal to the coherent detector 
by a few cycles per minute. ‘This consisted of a polyphase resistance capacity 
oscillator driving two balanced modulators in phase quadrature. The 20 kc/s 
signal was also fed to these two modulators in phase quadrature so that the sum of 
their outputs consisted of a single sideband differing from the original signal by 
the local oscillator frequency. This method was used in preference to a contin- 
uously rotating magslip because of the high frequency of the reference wave. 

As image frequency rejection was not applied, signals were received in two 
bands on either side of the local oscillator frequency, the total separation between 
these two bands being twice the intermediate frequency. ‘Thus, as the source 
moved through the aerial beams, two sets of fringes of slightly different frequency 
were produced. Hence if the source was observed for long enough these fringes 
moved in and out of phase, producing fringes with an amplitude modulation 
in the form of a rectified sine wave. Thus it can be seen that the rate of variation 
of the fringe amplitude with the delay was proportional to the separation of the 
two radio frequency passbands. This meant that, in order to keep the time 
delays as non-critical as possible the rather low intermediate frequency of 
13°5 Mc/s had to be used. Also to ensure that the phase dispersion of the 
two channels remained as nearly the same as possible most of the amplification 
was obtained in an amplifier common to both channels. This resulted in a 
sacrifice in signal to noise ratio because the signals from the two aerials had to be 
multiplied together in a detector instead of a linear multiplier. 

Operational method.—After the baseline had been changed, the movable 
aerial was levelled and the baseline surveyed with a theodolite, so that the variation 
of the difference in delay between the two signal paths could be calculated as a 
function of the hour angle for each source. This function was then used to 
tind the length of cable that had to be added to one of the intermediate frequency 
cable links to compensate for the delay. ‘The first derivative of this function 
divided by the wavelength gave the fringe frequency that would have been 
observed in the absence of the phase shifter, so that the apparatus could be set 
to give the required fringe frequency. This frequency was chosen to be low 
enough for the attenuation produced by the 12 sec time constant in the d.c. 
amplifier to be small and large enough for a number of complete cycles to be 
obtained during a single transit of the source through the aerial beams. 

As the aerials could not be made to track on a source, they were set for a given 
hour angle and the source allowed to drift through the beams. By this method 
it was found possible to make a number of observations near to meridian transit, 
usually at intervals of twenty minutes. 

Since the beamwidth of the larger aerial only amounted to one degree and 
as the delay had to be set to the nearest metre, the effect of varying the altitude and 
azimuth of the aerials and the length of the delay cable from the calculated optimum 
values was investigated experimentally for each baseline. 

As far as total power was concerned the apparatus consisted of two straight 
receivers with their outputs added together. Thus it was not possible to record 
the changes in total power due to the source passing through the aerial beams, 
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because of apparatus instability. This meant that the absolute fringe visibility 
could not be measured, and so fringe amplitudes were measured, corrections were 
applied for gain changes in the equipment and then the results were normalized 
to give unit fringe visibility at zero baseline. 


TABLE I 
Cassiopeia (23N5A) 
E. W. Baselines N. S. Baselines 
Length of Normalized Length of Normalized 
baseline in fringe baseline in fringe 
wavelengths amplitude wavelengths amplitude 
120 0°98 + 0°05 136 0°89 +0°05 
142 0°99 430 0°77 
400 0°63 726 0°27 
550 0°47 880 0°06 
700 0°32 1730 0:08 
894 <o'1o 
1290 0'13 
1750 <o'lo 
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Fic. 2.—Cassiopeia (23N5A) East-West. Abscissa : length of baseline in wavelengths. 
Ordinate : normalized fringe amplitude. 
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Fic. 3.—Cassiopeia (23NsA) North-South. Abscissa : length of baseline in wavelengths. 
Ordinate : normalized fringe amplitude. 


The Cassiopeia (23N5A) source.—Measurements of the fringe amplitude of 
the Cassiopeia source have been extended in the East-West direction from 
700 wavelengths up to 1750 wavelengths. These measurements together 
with those published previously are shown in Table I and plotted in Fig. 2. 
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These confirm the original conclusion that the diameter of the source at this 
frequency is not larger, and not more than 25 per cent smaller, than at metre 
wavelengths. They also show that the second maximum in the Fourier transform 
is only 10 to 15 per cent of the primary maximum, which is taken to indicate 
that the source is devoid of very large scale fine structure such as is found in the 
Cygnus source. 

The fringe visibility has also been measured with approximately North-South 
baselines of lengths up to 1750 wavelengths. The results of these measurements 
are shown in Fig. 3. It can be seen that the diameter here is the same as the 
East-West diameter to within the errors of measurement, thus showing that the 
source is nearly circular as it is at longer wavelengths. Here again there is no 
pronounced second maximum in the transform. 

The Cygnus (19N4A) source.—As it was suspected that the Cygnus source 
might be markedly asymmetrical with its minor axis in an approximately North— 
South direction, measurements were taken first with North-South baselines in 
order to find the direction of the minor axis accurately. This was done by 
observing the source for several hours on either side of upper culmination to 
measure the variation of fringe amplitude with hour angle. 

As the source was only 12° south of the zenith at transit and moved so as to 
describe a circle around the pole star, baseline foreshortening was a maximum at 
transit and so was never more than 3 per cent. Thus any variation of fringe 
visibility with hour angle was either instrumental or was due to variations in the 
angle @, between the source meridian and the baseline projected on to the celestial 
sphere at the source. 

Systematic variations of fringe amplitude with hour angle due to receiver 
sensitivity changes were not considered likely because each set of observations 
was taken within a few hours and the mean of a number of sets of observations 
obtained, while delay variations were compensated for by adding calculated 
lengths of cable to one channel for each individual transit of the source through 
the aerial beams. The attenuation produced by this cable was allowed for 
numerically. 

As the hour angle of the source varied, the elevation changed and, in view 
of the short wavelength, it was possible that the aerial gains might change. As 
total power from the source could not be recorded from either aerial this possible 
systematic error could not be eliminated directly. Thus measurements were 
made with baselines at different position angles so that a given value of # occurred 
for different hour angles from one baseline to the next. 


TaBLe II 


Minor axis of Cygnus (19N4A) 


Length of Normalized 
baseline in max. fringe 
wavelengths amplitude 
136 0°99 tO°IO 
430 1*00 
726 0°88 
1730 0°97 
2480 o'82 
3380 0°75 





No. 1, 1959 radio sources at 10°7 cm 31 


Figure 4 shows the fringe visibility plotted against hour angle and values of 
6 for three different baselines. It can be seen that the primary maximum occurs 
at a constant value of @ but at different hour angles for each baseline. Also it 
is evident that the fringe amplitude changes more rapidly as the length of the 
baseline is increased. As the maximum fringe amplitude occurs when the minor 
axis is in the same direction as the projected baseline, the position angle of this 
axis is 19 + 2°, so that of the major axis is 109° + 2°. 

It is interesting to note that the results for the longest baseline have a second 
maximum. ‘This can readily be explained in terms of a source model consisting 
of two centres of emission similar to that suggested by Jennison and Das Gupta 
(1953). The curves drawn in Figs. 4, 5 and 6 are the calculated fringe visibility 
for such a model. 
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Fic. 4.—Position angle of the minor axis of Cygnus (19N4A) Abscissae: top, 0 in degrees 
common to all diagrams ; others, hour angle in hours. Ordinates : normalized fringe amplitude. 

In Fig. 5 the maximum fringe visibility against length of baseline is plotted. 
This shows that the fringe visibility is only just beginning to fall at the longest 
baseline, from which it is estimated that the width of the minor axis is only about 
20”. 

In obtaining the Fourier transform of the intensity distribution along the 
major axis a difficulty was encountered in that the ground was not available for 
baselines of the necessary length in the required position angle. Thus measure- 
ments were made with baselines in position angles of approximately 80° and 
observations were in general taken several hours before transit. The effective 
length of baseline for each observation was then calculated by multiplying the 
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actual length by the sine of the angle between the aerial beams and the baseline 
to correct for foreshortening, and by the cosine of the angle between the major 
axis of the source and the baseline projected on to the celestial sphere to allow for 








© 400 800 1200 1600 2000 2400 2800 3200 300 4000 
Fic. 5.—Cygnus (19N4A) minor axis. Abscissa: length of baseline in wavelengths. 
Ordinate : normalized fringe amplitude. 
Tas_e III 
Mayor axis of Cygnus (19N4A) 
Actual length Effective length Normalized 
of baseline in of baseline in fringe 
wavelengths wavelengths amplitude 
180 180 o°96+0°'10 
890 630 0°57 
800 0°30 
1300 1020 <O'15 
1180 0°24 
1750 1500 0°68 


1670 0°67 


\ /. " 


° 20 400 600. 800 1000 1200 1400 1800 1800 2000 








Fic. 6.—Cygnus (19N4A) major axis. x Jennison and Das Gupta. © 10 cm results. 
Abscissa : length of baseline in wavelengths. Ordinate : normalized fringe amplitude. 


loss of resolution due to the projected baseline not being parallel to the major 
axis of the source. ‘This second correction was considered justified since it had 
already been shown that the lengths of the baselines used in the East-West direction 
would be insufficient to produce significant resolution in the direction of the 
minor axis. When these corrections were applied, it was found that the maximum 
resolution was obtained shortly after the angle between the projected baseline 
and the major axis of the source had passed through zero, since here the increasing 
angle offset the increasing resolution due to decrease of baseline foreshortening. 
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In Table I11 two values of fringe amplitude are given for each of the longer 
East-West baselines. One corresponds to the maximum resolution obtainable, 
while the other is the mean of the results taken on either side of the maximum. 
Figure 6 shows these results plotted against the effective length of baseline. For 
comparison the results given by Jennison and Das Gupta (1956) for baselines 
with position angle 113° are plotted on this diagram. It is seen that the first 
minimum in the Fourier transform occurs at appreciably shorter baselines and 
that the second maximum is noticeably larger at the higher frequency. If these 
results are interpreted as indicating a source with two separate centres of emission, 
this would mean that the two parts of the source, while remaining much the same 
size, had increased their angular separation by about 20 per cent. 

Discussion.—These results suggest that the size of the Cygnus source has 
increased in going from a frequency of 125 Mc/s to 2800 Mc/s while Conway 
(1956) comes to the conclusion that no appreciable change takes place in going 
from 125 Mc/s to 500 Mc/s. However, on examining his graph for the Cygnus 
source it is found that he has plotted the results of Jennison and Das Gupta as 
though they were fringe visibility whereas in fact they are equal to the square 
of the fringe visibility (Hanbury Brown and Twiss 1954). When the square 
root of Jennison and Das Gupta’s results are plotted together with those due to 
Conway it is again found that the Cygnus source appears larger at the higher 
frequency. 

Taking the major axis of the source to be in the East-West direction, Jennison 
and Das Gupta (1953) found reasonable agreement between their results and 
those of Smith. However, if these two sets of results are compared on the 


assumption that the major axis has position angle 109° at metre wavelengths, 
there is again an indication of increasing diameter with increasing frequency, 
although here the two frequencies differ by a factor of less than two. 
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VELOCITY ELLIPSOIDS AND THE GRAVITATIONAL 
POTENTIAL OF THE GALAXY 


P. A. Wayman 
(Received 1958 December 12) 


Summary 


In 1956 a model of the gravitational potential of the Galaxy was derived 
by M. Schmidt from an analysis of the circular velocities indicated by the 
21 cm radiation from interstellar material and of other data. It is found 
that, with a little modification, Schmidt’s model can be regarded as an 
example of a type of steady-state galaxy described in 1941 by Camm. The 
velocity ellipsoid at the Sun’s radius in this type of galaxy may have three un- 
equal axes and it is possible to discuss the variation of these axes in different 
parts of the Galaxy. One result is a tilt of the ‘* direction of star-streaming ”’ 
a fairly short distance above or below the galactic plane. 





1. Introduction.—Certain features of the Galaxy, such as the presence of 
dissolving clusters and the location of certain types of stars in spiral arms, deny 
the possibility of representing the Galaxy entirely as a system of stars in a steady 
state, and yet the age of certain other members indicates that the general shape 
and size of the Galaxy is unlikely to be changing rapidly. Oort (1) has thus en- 
visaged the formation of a disk-structure as representing a quasi-stationary stage 


in the condensation of the interstellar material towards the centre of the Galaxy 
wherein further collapse is offset by the galactic rotation; groups of interstellar 
material have interacted among themselves until they all partake in rotation 
about an axis which was originally the axis of the total angular momentum of the 
material concerned, a disk-structure having minimized this interaction. Con- 
version of this material into the older background stars of the galactic disk, 
following the formation and subsequent dissolution of clusters, may be expected 
to produce amongst these stars some features associated with a steady state in 
their distribution and their velocities even though the relaxation time of their 
own mutual interaction may itself be far too great to have caused such features. 
At the same time it may also be possible to trace the remnants of group motions 
among the same stars, as has recently been re-emphasized by Eggen (2). In 
any Case it is important to know what features of the Galaxy are, and what features 
are not, those which might be found in a steady-state galaxy. It is in this spirit 
that we consider here some features in the relation between the distribution of 
stellar density and the observed motions of stars in the region near the Sun, 
taking into account what is known of the gravitational potential of the Galaxy as 
a whole. 

For this purpose a suitable model of the distribution of mass in the Galaxy is 
one with rotational symmetry about an axis and with a plane of symmetry, the 
galactic plane, perpendicular to this axis. The representation of stellar motions 
by the ellipsoidal hypothesis of Schwarzschild (3) was originally empirical, but 
it was soon realized that it offered some chance of establishing possible relations 
between stellar motions and stellar density. Among the first attempts to do this 
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was that of Jeans (4) in 1915. One of Jeans’ results was that in a system with 
axial symmetry a non-isotropic distribution of velocities was only possible, for a 
steady state, if the distribution of velocities at a given point, or in a small region, 
had symmetry about an axis perpendicular to the plane through the point con- 
sidered and containing the axis of symmetry of the system as a whole. That is, 
star-streaming, if it exists, must be in a ‘‘transverse’’ direction. Reference to 
this result has frequently been made, up to the present time, (see, e.g., Smart (5), 
Oort (6), Schmidt (7), Woolley (8)), even though other possibilities have been 
discussed, first by Eddington (9), writing simultaneously with Jeans, and later 
by Clark (x0), Chandrasekhar (x1), Camm (12) and Parenago (13). It was 
curious that Eddington’s result predicted the existence of radial star-streaming 
for the case of axial symmetry and some time later Eddington himself (14) 
thought that between them he and Jeans had shown that neither radial nor 
transverse star-streaming were possible in the axially-symmetric steady-state 
galaxy. Chandrasekhar (15), however, showed that Eddington’s formulation, 
though valid, was not general; it had seemed too restrictive by the forms that it 
imposed on the gravitational potential function and Chandrasekhar pointed out 
that Eddington’s solution involved the unnecessary assumption that a function 
of three independent variables would be a perfect differential. 

In the meantime the study of motions within the galactic plane based upon 
the assumptions of a steady state had been shown by Oort (6) and Lindblad (16) 
to yield important results connected with the newly-discovered existence of 
galactic rotation about a distant centre. These results were: 

(i) a striking interpretation of Stromberg’s empirical parabolic relation (17) 
between velocity dispersion and mean motion for different groups of objects in 
the Galaxy ; 

(ii) an interpretation of the strong asymmetry existing for the velocities of 
the high-velocity stars; and 

(iii) a relation between Oort’s constants of galactic rotation, A and B, and 
the ratio of axes of the velocity ellipse in the galactic plane. 

Besides these important and well-known developments, Oort (18) had 
carried out an investigation of stellar motions in a direction perpendicular to the 
galactic plane based upon the assumption of a steady state existing in the density 
and velocity distributions in this direction for the region at the Sun’s distance 
from the galactic centre. 

It became generally accepted that Jeans’ result concerning the equality of 
the ‘‘radial’’ and ‘‘z’’ axes of the velocity ellipsoid was inevitable under the 
assumption of a steady state unless the potential function could be separated. 
That is, unless 6*/¢Réz=0. This condition would hold at the galactic plane 
but not over an extended region. A more complete picture of the consequences 
of a steady-state hypothesis was given by Camm (12) and it is of interest to know 
how far his results can be applied to the information on the density distribution 
in the Galaxy which has been provided by the interpretation of the 21 cm radiation 
profiles in terms of circular velocity w,(R). Kwee, Muller and Westerhout (19) 
published a discussion of such profiles and Schmidt (7) utilized these data 
to the full in constructing a potential function based on a model consisting 
of a series of spheroidal mass distributions. He also made use of information 
on the potential gradient in the z-direction at the Sun’s distance from the 
galactic centre. Schmidt’s model seems to be restricted by having te 


3* 
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assume uniform eccentricity inside his elemental spheroidal mass-distributions. 
This limitation shows itself in giving what is perhaps a too great degree of flat- 
tening for the central regions of the Galaxy and possibly leads to the postulation 
of a considerable amount of unidentified matter contributing to the mass of the 
Galaxy. Nevertheless Schmidt’s ‘Table 7 seems to provide the best available 
numerical values of the potential function if we bear in mind these limitations, 
as well as those which stem from the uncertainties in the basic data; for example, 
in the value assumed for the Sun’s distance from the galactic centre, 8-2 kpc. 

2. Steady-state theory in cylindrical coordinates.—If (q,, q2, q3) are generalized 
coordinates of position with canonical momentum coordinates (f,, po, p3) and 
there exists a frequency function for the distribution of stars in position and 
velocity 

AN ="E" (415 925 45» Pis Par Ps) 491 442 493 dP, dps aps 
Liouville’s theorem, or the equation of continuity, applied to the six-dimensional 
phase-space, states that ‘Y’ must satisfy the equation 


OF’  ~ foH ot’ cH ct’) 
—+Sis-a——s-a— mo (1) 
ot s (Op, 9q, 09,9, J 
where H is the Hamiltonian function per unit mass and the p’s are generalized 
momenta per unit mass. 

H=T (Px, qs) + P(4s) (2) 
where 7’ is the kinetic energy per unit mass and ® the negative potential energy 
per unit mass. 

It was shown quite generally by Camm (12) and Chandrasekhar (20) that 
for stellar systems which are in a steady state and which are of finite spatial extent 
it is only possible to have differential motions if the potential is characterized by 
axial symmetry. We consider this to be a suitable representation of the Galaxy 
and we choose cylindrical coordinates 


1 =R, q2=9, 93=3. 


The galactic centre (Fig. 1) is C, taken to be a distance R,=8-2kpe from the 
Sun in a direction which has galactic coordinates /= 3277-5, b=o°. 


@ 


Z 
4 


2 P 





Fic. 1. 


Axial symmetry in © implies 0@/0@=o0 and if, in addition, the plane z=o0 is a 
plane of symmetry, D=@(R*, z*). The kinetic energy per unit mass is 
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T= 4(Tl? + ©?+ Z?) where II, ©, Z are velocities at P (Fig. 1) in the “‘ radial’’, 
“‘transverse’’, and ‘‘s’’ directions respectively. In terms of q,, qa, 43; 

T=} (G1? + 91742? + 9s") 

oT ; 


Then p= 24, 


=q,=1 


oT 
p2= Oge = *g.= RO 


eS 
Ogs 
and H= (py +4 Pe +t) + (ay 4) 
If we introduce the observed sein distribution of velocities ‘’ where 
¥(R, 6, z, 1, 0, Z) R dR dé dz dil dO dZ 
=F" (Gis Ga» Ys» Pir Pas Ps) EP, Az dps dq, dq, dqs 
OY _ OY" | OY" apy _ 2" |g a" 
OR 04, op, OR dy ops 
_ OF" dp, _ pt” 


P3= jy =Z 





we have 


and 


ee ‘ a Y or 
Liouville’s theorem for ‘’ under the steady-state condition =0, = =O 


ot CT) 
becomes 
I ov ee on 2) or OOF a® oF _ 
OR dz oR R/ dal R00 dz OZ 
Under the generalized ellipsoidal hypothesis of Chandrasekhar we postulate 
=‘'(O) where O=all? +50? — 2600, + cZ? + 2(f0Z + g11Z+Al10) +0. 

The coefficients a, b, Qo, c, f,g, hand o are functions of position only. We can write 
d¥ = =% dO and O must satisfy (4) for all values of II, @, and Z. The 
equation obtained by substituting O in (4) in terms of a, 6, etc., gives a polynomial 
of the third degree in II, @, and Z. The coefficients of this polynomial can then 
be separately equated to zero, if this equation is an identity for all values of 
Il,©@,Z. A number of equations that the coefficients a, 5, etc., in Q must 
satisfy are obtained and we can take over Camm’s result (12) for the solution of 
these equations for the present case. We must have :— 

A = Oy + a2” f=o 

b= % + a2" + aR? g= —a,Rz (5) 

c= Xs, + aR? h =O | 

= aR _ aR 
5b — (&o + a2" + aR?) 
do 


dr "oR 


da om = | 





(6) 


(7) 


, “2s 


where a, %1, %, %3, %4 are constants for the system of stars being considered. 
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From (7), equating 0?c/0Rdz in the two equations, it is necessary that ® 
should satisfy the equation 


( o® > ee eo ap a 
2543 5R — SRS + (3 RE) Se + Re ( )}+ (0-2) 


“oR Os oRos 


(8) 


OR? dz® 


This equation is equation 23 of Oort (6). He inserted the solution 
Xz =O, Aq = My 
presumably on the basis that this differential equation is too restrictive on ® 
otherwise. Camm showed, however, that if 
t,4#O and way, 
a solution exists of the form 
= (F(E)— G(m)i/(E—1) (9) 


where (€, 7) are elliptic-hyperbolic coordinates related to R and z by the equations 


R’=(&+B)(n+B)/28 and 2?=—(€—8)(n—8)/2B 
the constant 8 being B = (a) — «,)/2%5. 
Further, if y= }(a+2,), 
(da + yD} /a3 = {EG(q) — nF(E)}/(E—1) (10) 
F(€) and G(») are arbitrary single-valued functions of € and 7. 
For a group of stars governed by ’ as a frequency-function of the form 


Y =e-@ the stellar density v= { | | Wall dO dZ is 


11, 0, Z, 
v= 732(abc)-H2 @f-2+bOn), (11) 

If we consider only the points (R, 0) we can define the circular velocity 
©,.=Rw,, given by 0,2/R=(e@/eR), ,. At these points g=o0, so that, using the 
first of equations (7), 0,2=(R/2a)da/0R. Also, a=%, b=a49+%.R?, c= 44+ %,R? 
and b©,? = «,?R?/b so that 

2 (0,2) = 20,229) 
aR | 0°) = 20)7a5/R. 

Differentiating (11) w.r.t. R and inserting these expressions for a, 6, c, and 
0/2R(6©,?) we easily obtain 

r0y oR _ _ GR, M%1g2_ 90). (12) 
vOR a+aR% agt+a,R? R 
With a suitable change of notation this equation reduces to equation 31 of Oort 
(6). 

3. Schmidt’s model of the Galaxy.—The galactic model of Schmidt (7) is 
based upon the superposition of a number of non-homogeneous ellipsoids whose 
densities, varying in a particular way, are adjusted to fit primarily the circular 
velocities ©,(R) revealed in the inner parts of the Galaxy by the 21 cm radiation 
from interstellar hydrogen. For the outer parts an extrapolation was necessary. 
Existing information concerning the distribution of matter perpendicular to the 
galactic plane is also used. Schmidt finds that the superposition of four such 
non-homogeneous oblate spheroids gives a good approximation. The final 
solution employs the addition of a few homogeneous spheroids which are of 
quite small total mass. The final result is tabulated as —®, 0@/dR and 90/dz 
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in Schmidt’s Tables 7, 8, and g respectively. These data are used as the basis 
of the present investigation. Schmidt himself points out that his model suffers 
to a certain extent in having to assume uniform eccentricity inside his elemental 
spheroids, in that the resultant model is possibly too flattened towards the centre 
of the galactic system. However, lacking further information on the distant 
distribution of matter perpendicular to the galactic plane, no improvement 
seems possible. It is of interest in this respect that Parenago (13) suggested a 
form of ® given by ®=@(R).¢(z); he represented observed densities of 
RR Lyrae stars on the basis that (0, z)/@(o, o)= (Rp, z)/M(Ro, 0). There is 
however some uncertainty in his interpretation of the observations and Parenago’s 
form of ® is not in any way compatible with Schmidt’s model. It implies that 


(5), , oc D(R, 0) 


Oz? 


and this cannot hold everywhere as it results in much too large values of K, = d@/dz 
a short distance above and below the galactic plane at large distances from the 
centre. 

Schmidt attempts to identify some portions of his derived density distri- 
bution with observed groups of galactic objects, but there remains a considerable 
amount of mass to be ascribed to ‘‘ unknown objects’’. It is possible that 
errors in the adopted values of R, (8-2kpc) and the Oort constants of galactic 
rotation, A(+19°5km/seckpc) and B(—6-9km/seckpc) are connected with this 
anomaly. We shall here assume that Schmidt’s potential function is the best 
available and we shall endeavour to discover whether Camm ’s form (g) is 
suitable for the representation of Schmidt’s values. 

Kuzmin (21) obtains a solution to equation (8) which is identical to Camm’s 
if we transform (&, 7) to Kuzmin’s (€,, €) by the equations 


Z2= — 2B; fer 5 (E+); ét—1= 5 (n+). 


His evaluation of this solution as a galactic model is not, however, based on a 
complete ‘‘observational’’ model such as Schmidt’s but is more hypothetical. 
Further comparison of the present account with Kuzmin’s will be made later. 

4. An appropriate value of 8.—For points with R=o the coordinates (€, 7) 
are £= — 8, n=2?+8 and for points with z=o they are £=R*—f, n=8. It is 
easy to show that if D={F(£)—G(n)}/(€—7n) the following equation must be 
satisfied at all points (R, 0): 


(se). ~ =a" (aR), * B=aB(B*). 
+ pe g(0)-O(R,o)}}. (23) 


If we wish to apply this equation for R= R,=8-2kpce we can take the following 
values for the various differential coefficients from the data of Schmidt’s Table 7. 
Units are 100 km?/sec? or 100 km?/sec? kpc?. 
M(0, 0) = — 1436 ®(Ro, 0) = — 410 

ap 


Bree MGR (SB) 
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Inserting these values into (13) gives a quadratic for 8 with one root 8 = — 2-25, 
and a positive root which is rejected as it would involve difficulties at the point 
z=0, R*=28. 
A better value of 8 near to — 2-25 was then found in the following way. We 
have 
®(R, 0) ={ F(R’ — 8) — G(8)}/(R?— 28) | 
and (Ry, 2) =(F(E) — G(n)}/(E—2) 
where £+n=R,? +27 =67:24 +27 ( (14) 
and én = B(67:24-2*— 8). } 
Probably the best determined values of Schmidt’s Table 7 for ® are those of 
®(R, 0) for all R from o to 16-4kpe and of ®(R,, z) for z up to 18kpe. Using 
equations (14) and adopting trial values of 8 near to — 2-25, F(£) and G(n) could 
be deduced, thus giving 
F(—8)-G(2?+8) 
(0, z)= . 
(02) = +28) 
A comparison of these quantities could then be made with Schmidt’s values and, 
as shown in Fig. 2, the value B= —2-5 gave the best fit. Because the com- 
putation of G(y) depended upon the difference of two products, the solution 





log $ (0, z) 














06 0.7 
log (1 + 2) 
z inkpe 


Fic. 2. 


became insensitive for points (0, z) beyond s=1kpe and could not be expected 
in principle to give O(o, z) beyond z?= —2§8, or z=2-2kpc. The point R=o, 
z= 4/ —28 is a singularity in the (£, ») coordinates. 

Slight modifications were needed in Schmidt’s values of ®( Rg, z) for z > /4 kpc 
to give, with B= — 2-5, Schmidt’s values of (0, 2) over the range z=1 to 2kpc. 
A comparison of the two sets of ®(R,, z) is shown in Table I, while Table IT 
gives an abbreviated formof Schmidt’s Table 7 and its counterpart, using 
Camm’s form of ® evaluated as above with 8 = — 2°5. 

In Table II it can be seen that the region o< R<4kpc, z>0 is not satisfactorily 
represented. Inthe region of the singularity at (o, 2-2) some revision of Schmidt’s 
potential function is required before it will fit Camm’s form well, even though 





the gravitational potential of the Gaiaxy 


TaBLe I 


P(R, z) 
Schmidt 


100 km?/sec? 


(Ro, 2) 
modified 


100 km?/sec® 


Schmidt 
R (kpc) ° 
z (kpc) 
° 1436 
I 989 
4 539 


10 275 
18 160 


354 
332 
312 
273 
228 
180 
150 


TABLE II 


2°05 4°10 6°15 8:20 10°25 12°30 


1012 
851 
520 
270 
160 


748 
675 
460 
260 
158 


Camm Function with B= —2°5 


R (kpc) 
2 (kpc) 
° 1436 

I 980 

4 539 
10 275 
18 160 


° 


2°05 4°10 


1012 748 
764 674 
563 501 
274 256 
153 146 


554 
518 
395 
247 
154 


410 
396 
332 
228 
150 


354 
332 
310 
269 
227 
171 
136 


— O(R, z) Units: 100 km*/sec* 


318 
312 


258 
257 
280 239 
208 =: 188 
450 == 


12°30 16°40 


258 
254 
231 
179 
126 


188 
187 
174 
148 
116 
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elsewhere the fit is good. Kuzmin, by using a particular smooth function for 
density and potential, avoids difficulty with this singularity. 

Another check on the value of 8, this time for points (R, 0), may be made by 
considering how well equation (13) is satisfied. A comparison of the logarithmic 
values of 6°@/dz* found from (13) (filled circles) and from Schmidt’s Table 
direct (continuous curve) is shown in Fig. 3. The agreement is fairly good and 
neither here nor, for the most part, in Table IJ are the discordances considered to 
be large enough to say that the gravitational field of the Galaxy is not such that a 
system of stars can maintain quite closely a steady state in the way envisaged by 
Camm. 

5. The axes of the velocity ellipsoid—Camm pointed out that in the present 
analysis knowledge of ® does not imply the values of the axes of the velocity 
ellipsoid, which are still left with at least two degrees of freedom. It concerns 
rather their variation from point to point. Suppose that we have a set of stars 
in a steady state with axes for the velocity ellipsoid in the Sun’s neighbourhood 
(R=8-2kpc, z=0) 

o, = 41 km/sec 


o), = 25 km/sec 


? 


a, = 18 km/sec. 
These figures are taken from Allen (22), p. 221. a=1/20,”, etc., and, using 
equations (5) and B= }(%)—«,)/x3= —2°5, we obtain 
= 2° 1074 
1076 
5 so 


2 
3°83 10~% 


the units for x) and «, being sec?/km? and for x, and «3, sec?/km*kpc?. ‘These 
values (15) are satisfactory in that none are negative and the value of «, in (6) 
may be found by making the maximum value of ©, at points (R, 0) equal to the 
maximum value of ©. According to (6), for points (R, 0), 0O)/¢R=o where 
R=V %/2?=6-3kpe, from (15). ©, reaches a maximum quite near this value 
of Rand, with «, = 0°0212, Table III gives a comparison of ©, and 0... 


Tase III 

R ©, Oo 
kpc km/sec km/sec 
2°05 178°3 132°1 
3°075 198°4 179°4 
4°10 211°4 205°3 
6°15 224°6 224°6 
8-20 216°2 217°0 
10°25 192°1 200°3 
12°30 171°! 182°3 
16°40 142°0 150°3 


Assuming that these values of ©, are accurate, the consequences of the differences 
in Table III on a completely steady state picture are shown in Table IV, 
column 2. Using equation (12) the required values of (1/v) ov/OR are calculated. 
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If the system of stars with this velocity function are to give density gradients 
in accordance, roughly, with the values of Schmidt’s Fig. 1 we expect 


I Ov 
= oneee we = OG 
v 


over quite a large range of R, since the total density has to decrease by a factor e 
in about two parsecs; Schmidt discusses the rather high values of 1/v dv/OR that 


TasLe IV 
R I Ov, I Ov, 1 0M I Ov 
OR %.0R MéR v3 OR 
kpe 
2°05 —4°'28 — 5°34 —0°32 — 4°60 
4°10 —0°55 —I°I9 «=~ 0°32 —0'87 
6°15 —O'ls —0'50 — 0°32 —0°47 
8:20 —o'l4 —0°32 — 0°30 —0'44 
10°25 +0°035 —0'06 —o'18 —O'rs 
12°30 +0°056 0°00 —0°075 — 0°02 
16°40 —0°'021 — 0°04 —o'oIo — 0°03 
occur in his model. We can alter the assumed value of «, to remove the positive 
values in col. 2 of Table IV and obtain those values shown in col. 3. These values 
are not altogether satisfactory, especially for R<5kpc, and an attempt was 
therefore made to see whether a minor modification in the form of Y would give 
any chance of improving the situation. 


Suppose we take (II, 0, Z) =(V?—Il? — 0? — Z?) e-® where 
O=all* + 60? — 2b00, + 2gI1Z + cZ* +a, 


as before. We may identify V with the escape velocity given by V?= —2® so 
that ‘’=o when I]?+6?+Z?2=V?. It proved possible to maintain equations 
(5) and (6) as the steady-state solutions for a, b, c, g and ©, as the substitution of 
O’ = — In (V?—II*— ©? — Z?) + EC into (4) led to terms being added to the poly- 
nomial in II], © and Z (not given here) which were identically zero. An 
expression for the density was evaluated by making the approximation 


dll dO dZ= | | | Vall dO dZ. 
Il’?+@©?+ Zt=<vy 
The expression for v may be written, for points (R, 0), 


y= e(-2+O% M/V abc (16) 


where M=V7z(V2-0,2]— ar 5+} +3). 


In the interpretation of the observed velocities the nary of the constants 
a, b, c, g, and @, are the same as the corresponding quantities in the true ellip- 
soidal hypothesis, or very nearly so, since for small values of Il, ®—@,, and Z 
the quantity (V?—Il?— ©?— Z?) varies only slowly, even if V is no more than 
1°4@,). The fourth column of Table IV gives the values of 1/M 0M/dR and the 
fifth column the values of'1/v v/0R from (16). There remain the large negative 
values for R<5kpc and rather small values for R>1okpe but it is seen that a 
relatively small change in the nature of ‘Y may affect the satisfaction of this 
steady-state condition. 
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For other values of A= %) + 4,2" 

b = a + a2? + aR? 
C=a,+ agh? 
we may briefly consider the consequences that follow from the relation 
Ay — $23 = Hy 

together with the probability that «,/x%9~0-o3 (in the units here employed). 
The latter condition follows from the necessity for @, to reach a maximum some- 
where near R=6kpc. 

To maintain values of 1/v dv/dz of the order of — 4 along the line of points 
(o, z) it is necessary that there should exist, towards the central regions, a 
majority of objects with a,~4:10-°, or o,~110km/sec. We cannot from this 
put definite values to x and «, except as regards their probable order of magnitude. 
Using the units here employed we expect, numerically, from the form of a, 3, 
and c for large values of R and z, that o<a,~a3,<a, and hence «4—% = 5%, >0. 
We then deduce that («, + R,?«,) should be rather greater than «, but only slightly 
greater than («+ R,?x».) or, in terms of the velocity ellipsoid axes, o,; >og~¢z.- 
The solution x; =0, %=a,is nevertheless a possible one, giving o,;=0,>0@, as 
in Jeans’ solution, though of course the analysis of Camm’s solution then does 
not hold and no value can be assigned to f. 

The distribution of the velocities of high-velocity objects is so little known 
that it does not seem profitable to follow the various possibilities further. 

Taking the four values (15) of %, %, %; and a, the velocity ellipse in the 
plane 6=constant at a point R= Ry, z=0°5 kpc is 

3701311? — 1-41411Z + 15-432? = 10%. 

The non-zero value of the coefficient of II1Z in this equation gives a tilt of the 
axis 3° 15’ away from the line parallel to the galactic plane. In order to detect 
such a small tilt it seems that one would need the radial velocities of about 1000 
stars at approximately this distance whose proper motions are sufficiently well 
known for them to be assigned as having a positive or negative [l-value with 
some certainty. This is rather far beyond the resources of present catalogues. 

6. Conclusions.—The situation seems to be as follows. Schmidt has built up 
a model of the Galaxy to fit available data, assumed to be correct. He has 
constructed it in a semi-empirical way out of relatively few non-homogeneous 
spheroids of a particular kind of mass-law, superimposing a number of homo- 
geneous spheroids with a small total mass (1 per cent of whole). Camm’s 
functional form for a steady-state galaxy cannot be made to fit analytically the 
potential function within these spheroids but it is suggested here that the final 
numerical values for the potential function arrived at by Schmidt do fall roughly 
into Camm’s form, at any rate over the majority of the volume occupied by the 
Galaxy. ‘Thus it is suggested that the Galaxy has reached a shape in which a 
steady-state distribution of stellar density and velocity is possible to a fairly high 
degree of closeness and it has been shown that the ratio of axes of the velocity 
ellipse for stars in the solar neighbourhood is such as would be expected for such 
a distribution. Kuzmin (23) arrives at a similar conclusion by a slightly different 
process. His idealized models are effectively such as to allow a steady-state 


distribution but it is not clear whether his superposition of spheroids does in fact 
allow it. 
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The difficulty of satisfying Poisson’s equation, discussed analytically by 
Camm, has not, so far, been mentioned. It is impossible for a group of stars in 
a steady state with a uniform mass-distribution to produce a gravitational field 
in which a steady-state motion of the type discussed here is possible. The only 
explanation, as applied to the Galaxy, is to assume that the superposition of 
different systems within the Galaxy with different mass frequency-distributions 
nevertheless does allow a steady state. Given what is known of the potential 
function of the Galaxy, non-equality of the radial and z axes of the velocity 
ellipsoid near the Sun is not inconsistent with the steady-state hypothesis. 
Other observed features, of course, are, but they are more concerned with details 
than with the overall gravitational field. It is important in investigating non- 
steady features that we should be clear as to what features do in fact represent 
departures from the steady-state hypothesis. 
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THE RADIO EMISSION FROM THE DIRECTION 
OF THE ‘‘SUPERGALAXY’”’ 


Jj. R. Shakeshaft and F. E. Baldwin 


(Received 1958 December 29) 


Summary 


Observations have been made with the Cambridge radio telescope at a 
wave-length of 1°9m from which radio isophotes of the region of sky 
10°14" have been derived in order to investigate the features which have 
been ascribed to radiation from the “ local supergalaxy’’. The distribution of 
emission is not closely correlated with counts of galaxies. The observed 
intensities are several times greater than those to be expected from normal 
galaxies and counts of radio sources show that this factor cannot be explained 
by a small proportion of the galaxies having abnormally high radio luminosities. 
The radio brightness of one of the features is very large compared with the 
total background of extra-galactic emission. It is suggested that the radio 
features are not in fact associated with the “ supergalaxy ”’. 





1. Introduction —Surveys of the background radiation by Hanbury Brown 
and Hazard (1953) at 159 Mc/s and Kraus and Ko (1953) at 242 Mc/s have 
revealed that there is a bright band of radio emission running roughly 


perpendicular to the galactic plane at about 12" R.A. These authors have 
suggested that this band represents the integrated emission from the concentration 
of bright galaxies lying along a great circle crossing the galactic plane at longitudes 
l=105° and 285°. 

De Vaucouleurs (1953, 1956, 1958) has presented evidence to show that 
these galaxies and our own are members of a flattened system or ‘‘cluster of 
clusters’’ with its centre near the Virgo cluster. This he terms the Local 
Supergalaxy, after Shapley. In his later paper the radio data are summarized 
as additional support for this picture, but it is the purpose of the present paper 
to show that there are certain difficulties in such an interpretation of the radio 
feature. 

2. The shape and position of the radio feature in relation to the counts of 
galaxies.—One of the four sections of the Cambridge four-aerial interferometer 
(Ryle and Hewish 1955) has recently been used for a survey of the background 
radiation at a frequency of 159 Mc/s. The beam-width between half power 
points is 1°+2 in right ascension and 7° in declination. Contours of brightness 
temperature derived from these measurements are shown in Fig. 1, plotted 
together with the positions of galaxies from the Shapley-Ames catalogue (1932). 
This catalogue contains the positions of galaxies brighter than about mpg= 13 
and is believed to be substantially complete to that level. 

In the map of Fig. 1 the areas of greatest brightest temperature are shown 
hatched. The contours are noticeably affected by the presence of two bright 
radio sources. One of these, at «12"28™, §+12°, is Virgo A and the other is 
Cassiopeia A which, observed at lower culmination in a 1 per cent sidelobe of 
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47 
the declination reception pattern of the aerial, causes the long ridge extending 
from «12. 10™, 8+ 50° to ar1®50™, §+08°, 
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Fic. 1.—Contours of aerial temperature at 159 Mc/s plotted together with the positions of 
galaxies from the Shapley—Ames catalogue. |The contour interval is approximately 8°K. 


Apart from these sources, the main features to be seen on the map are: 


(i) A general increase in brightness from « 10" 30™-14" 00™ associated with 
the extended component of the galactic radio emission. 


(ii) An intense ridge between 5+08° and +27° extending from « 14"00™ 
to 12° 30m, 


This is part of a much larger feature, detected in many surveys, 
which runs from the galactic plane near /=o° towards the north galactic pole. 
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(iii) A bright region between 6 —15° and +05° and « 12"15™ to 13" 00", 
It is this feature that Kraus and Ko (1953) believe to be associated with emission 
from the Supergalaxy. 

(iv) A comparatively faint ridge of emission running roughly north and 
south between « 11"00™ and 12"30™ at declinations north of 6 +15°. Part 
of this feature is that studied by Hanbury Brown and Hazard (1953). 

Since faint features are not easily visible on a contour map if there are 
appreciable gradients in the background, sections across the map at several 
different declinations are plotted in Fig. 2. 

The enhanced radio emission in the regions (iii) and (iv) above can be seen 
to coincide roughly with the distribution of the bright galaxies, but the detailed 
correlation between the intensity of the radio emission and the density of the 
galaxies is poor. In particular the following points may be noted: 

(a) he area most densely populated with galaxies, centred on about 
z 12" 30™ and 6 12’, in which the Virgo cluster is situated, does not correspond 
with a peak in the radio isophotes. 

(6) The pronounced maximum of the radio isophotes at x 12" 40™, 6 —08 
does not seem to coincide with a high concentration of galaxies. De Vaucouleurs 
(1956) shows that obscuration can have little effect on the apparent distribution 
of the bright galaxies in this region of the sky. 

(c) The latter maximum seems, furthermore, to be connected to the ridge 
of high intensity (11) which shows no correlation with counts of galaxies. There 
is strong evidence for associating the latter feature with the Galaxy (Baldwin 
1955). 

The galaxies plotted in Fig. 1 are mainly those of the Local Supergalaxy 
having a high intrinsic luminosity and an apparent magnitude less than +13. 
It would be interesting to see if the radio isophotes were in better agreement 
with the distribution of fainter galaxies. There is an excess in numbers of 
galaxies in the north galactic hemisphere over the south, persisting down to 
an apparent magnitude of +18, which de Vaucouleurs attributes to members 
of the Local Supergalaxy with low intrinsic luminosities. Unfortunately, little 
work has been published on the distribution of galaxies in different magnitude 
ranges for the region covered in Fig. 1 and such information is available only 
for a zone from 6 +41° to +46. Data for this zone are given by Shapley and 
Jones (1940) in the form of the number of galaxies in nine square degrees of sky 
centred on various right ascensions along declination + 43°30’. ‘These numbers 
for three different magnitude ranges are plotted in Fig. 3 together with the aerial 
temperature measured at a wave-length of 1-9 m with a reception pattern centred 
on declination + 42°14’. Most of the variation of the radio intensity is caused 
by the main component of the galactic emission, but the feature under consideration 
can be seen at « 12", It will be noticed that the belt of excess radiation coincides 
with one of the smaller increases in number density of the fainter galaxies. ‘The 
emission associated with the considerably larger peaks in the galactic counts at 
13" and 16" is certainly not greater than 20 per cent of that at 12". 

Since the aerial beam has a width to half-power points of 7° in declination 
it extends beyond the limits of the areas of nine square degrees used for the 
galaxy counts. The counts of galaxies in areas of 25 square degrees given by 
Shapley and Jones were therefore compared with the radio observations. The 
results were substantially the same as those shown in Fig. 3. 
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Fic. 2.—Variation of the aerial temperature with right ascension for a number of different declinations. 
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3. The radio emission expected from the Supergalaxy.—Hanbury Brown and 
Hazard (1953) have calculated the intensity of the integrated emission which 
would be expected from the normal galaxies in a region near «= 12", = + 43°. 
They used the data given by Shapley and Jones, assuming that all galaxies have 
the same ratio of radio flux to light flux as that found for spirals not emitting 
enhanced radiation. 
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Fic. 3.—This figure shows the numbers of galaxies in g square degrees to different limiting 
magnitudes in areas of sky centred on 8+ 43° 30° at various right ascensions. Also plotted is the 
aerial temperature at 8+42° 14’ as a function of right ascension. 


The computed value of the integrated emission was 1-7 magnitudes (4:8 x ) 
less than that actually observed. 

Observations by Mills (1955) indicated that the emission from elliptical 
galaxies is considerably less than that from spirals and this result therefore 
increases the discrepancy. The most recent work on the flux densities of normal 
galaxies presented by Hanbury Brown and by Mills at the Paris symposium on 
radio astronomy, 1958, suggests that an average value for the radio magnitude 
at 160 Mc/s, m,(160) of a spiral galaxy with photographic magnitude mpg may 
be taken as m,(160)—mpg>0o where m, is defined by 

m,(160)= — §3°4—2°5 logy, S(160). 

Using this figure to derive the expected brightness temperature at 160 Mc/s 
from normal galaxies down to mpg 17°5 in the region near a12", 6 +43° a 
value of 5 °K is obtained, i.e. less than 20 per cent of that observed. Because 
some proportion of the galaxies will be ellipticals this estimate should be an 
upper limit. 

Although the observed intensity cannot be accounted for by the integrated 
emission from normal galaxies, there is the possibility that the radiation arises 
in abnormal galaxies having a much higher ratio of radio flux to light flux. For 
instance, since the collision rate of colliding galaxies depends on the square of 
their space density, a larger proportion of colliding galaxies might be expected 
in the Supergalaxy than in neighbouring areas of the sky. 

In this case, if the observed excess brightness temperature in the belt is to 
be explained as the integrated emission from such sources, then a consideration 
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of the actual numbers of radio sources observed in the region enables a lower 
limit to be placed on their space density. 

Suppose that the maximum excess brightness temperature of 80°K is 
the integrated effect of discrete sources each having an equal luminosity 
Pw.ster~1(c/s)~! with a space density of pMpc~*, then the flux density per 
steradian from distances between R, and R, is 


R, 2p 
| pt dr=pP(Ry—R,) 
R, r 


and the brightness temperature corresponding to this is 
2 
T= 5 pP(R.— Ri). 

Taking (R,—R,) as <15 Mpc, the approximate extent of the Supergalaxy 

according to de Vaucouleurs 
Pp>1-4x 10~ MKS units. 

The survey of radio sources at a frequency of 159 Mc/s recently completed 
at Cambridge shows five sources having a flux density > 15 x 10-*6 w.m.~* (c/s)~! 
in the area of sky 12" 00™<a<13"30™, —10°<5<+4° (approximately 
i steradian). If all of these were sources in the Supergalaxy* then }pR*=5 x 12 
where R is related to P by P/R?=15 x 10-*6 w.m.~? (¢/s) 7}. 

Thus ’ 

pP3? = 10% MKS units 
whence 
p> 700 Mpc-%, 
which is much greater than the density of ordinary galaxies, so that the radio 
isophotes cannot be explained in terms of galaxies having abnormally strong 
radio emission. 

This calculation ignores the effect of a dispersion in absolute luminosities 
which would lead to a still greater density of abnormal sources. 

Shklovsky (1954), using different arguments, also came to the conclusion 
that emission from individual sources could not account for the feature; he 
suggested that the emission was due to relativistic electrons accelerated by 
inter-galactic magnetic fields of the order of 10~* gauss throughout the whole 
volume of the Supergalaxy. This is not a hypothesis that can be easily tested 
but, in common with all the other theories that the radiation is associated with 
the Supergalaxy, it is, as will be shown later, difficult to reconcile with the total 
amount of the extra-galactic radio emission. 

Some further experimental evidence is available which has a bearing on this 
matter. De Vaucouleurs (1956) shows that there are other ‘‘supergalaxies ”’ 
not far from the local one. The nearest of these appears as an elliptical swarm 
of galaxies running from « 3"30™, 6 —7° to a 5", 6 —60°. He suggests that 
it is seen edge-on and that it is perhaps rather smaller than the Local Supergalaxy. 
Bolton and Westfold (1950) have surveyed this region at a frequency of 100 Mc/s 
with a 17° beamwidth and an estimated accuracy of +20°K. The swarm lies 
right across a minimum of their isophotes but no irregularities can be seen. 
The angular dimensions of the swarm are 10° x 60° and so an aerial temperature 
of about 100°-200° K might be expected. 

* Since the Supergalaxy is not apparent on source maps only 1 or 2 of the 5 could be associated 
with it. 


4* 
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The central part of the swarm is, unfortunately, not visible from Cambridge, 
but some of the scans at a frequency of 159 Mc/s cross it at a place where the 
concentration of bright galaxies is still marked. Fig. 4 shows the variation of 
aerial temperature with right ascension at 6 —28°. The position of the super- 
galaxy in question is also indicated and it can be seen that the upper limit to the 
aerial temperature from it is about 5 °K. 
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Fic. 4.—The variation of temperature at 159 Mc/s and 8—28° with right ascension. Also 
marked is the position of a swarm of galaxies noted by de Vaucouleurs. 


An examination may also be made of dense clusters of galaxies which, on 
Shklovsky’s hypothesis, should be very intense sources. The most readily 
detectable should be those in Perseus and Coma. Data on the numbers and 
types, e.g. spirals, ellipticals, etc., of galaxies in these clusters are meagre and 
it is not possible to calculate accurately the flux expected if there are no sources 
with enhanced emission. For a very rough estimate let us assume that there are 
1000 galaxies in each cluster, with a mean mpg of 16. Assuming also that the 
relationship m,,(160)—m,,=o may be applied, then the expected flux density 
at 160 Mc/s is 16 x 10°" w.m.~* (c/s)~}. 

If this value is compared with the observed flux density at 81-5 Mc/s for the 
Perseus cluster of 45 x 10-76 w.m.~? (c/s) ~! after the contribution from NGC 1275 
has been subtracted (Baldwin and Elsmore 1954), and of 10-7 x 10-** w.m.~?(c¢/s)~! 
at 400 Mc/s for the Coma cluster (Seeger, Westerhout and Conway 1957), the 
discrepancies do not seem to justify the belief that there is a large excess of 
radiation from these clusters. 

4. The integrated extra-galactic component.—In the brightest part of the 
radio feature the excess radiation (about 80°K at 160 Mc/s) amounts to more 
than 20 per cent of the total intensity in that direction. Most of the background 
emission is, however, galactic in origin. An upper limit on the amount of the 
extra-galactic component has been set by Baldwin (1955) from an analysis of the 
angular distribution of the galactic background. This limit is 500 °K at 81-5 Mc/s. 
Although the spectrum of the extra-galactic component is at present unknown, 
we can use the 81-5 Mc/s observations to show that the excess radiation from the 
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direction of the Supergalaxy is of the same order of magnitude as the total 
extra-galactic radiation. 

If the radiation originates in the Local Supergalaxy, then the first 15 Mpc 
of space measured in a line of sight from the Earth in the direction of the 
Supergalaxy must be intrinsically different from, say, the next 1000 Mpc. 
In view of the large scale uniformity found among the more distant galaxies, 
such a state of affairs seems unlikely. 

5. Conclusions.—The evidence presented in this paper does not support the 
physical association of the radio feature with the Local Supergalaxy and it seems 
possible that it originates within the galaxy. It lies close to the much more 
intense belt of radiation which is very probably galactic. The radiation from 
this latter feature has been discussed by Mrs H. Tunmer (1958). 

Our thanks are due to Mr D. O. Edge for allowing us to use information of 
the 3C survey of radio sources prior to publication. ‘This work was supported 
by a grant from the Department of Scientific and Industrial Research. 
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The central part of the swarm is, unfortunately, not visible from Cambridge, 
but some of the scans at a frequency of 159 Mc/s cross it at a place where the 
concentration of bright galaxies is still marked. Fig. 4 shows the variation of 
aerial temperature with right ascension at 6 —28°. The position of the super- 
galaxy in question is also indicated and it can be seen that the upper limit to the 
aerial temperature from it is about 5 °K. 


Extent of 





“Swarm of Goloxies” 


Brightness Temperature (Deg. K) 








0400 0300 
Right Ascension 


Fic. 4.—The variation of temperature at 159 Mc/s and &—28° with right ascension. Also 
marked is the position of a swarm of galaxies noted by de Vaucouleurs. 


An examination may also be made of dense clusters of galaxies which, on 
Shklovsky’s hypothesis, should be very intense sources. ‘The most readily 
detectable should be those in Perseus and Coma. Data on the numbers and 
types, e.g. spirals, ellipticals, etc., of galaxies in these clusters are meagre and 
it is not possible to calculate accurately the flux expected if there are no sources 
with enhanced emission. For a very rough estimate let us assume that there are 
1000 galaxies in each cluster, with a mean mpg of 16. Assuming also that the 
relationship m,(160)—m,,,=o may be applied, then the expected flux density 
at 160 Mc/s is 16 x 10°78 w.m.~* (c/s)~1. 

If this value is compared with the observed flux density at 81-5 Mc/s for the 
Perseus cluster of 45 x 10-78 w.m.~? (c/s) ~! after the contribution from NGC 1275 
has been subtracted (Baldwin and Elsmore 1954), and of 10°7 x 107° w.m.~*(¢/s)~! 
at 400 Mc/s for the Coma cluster (Seeger, Westerhout and Conway 1957), the 
discrepancies do not seem to justify the belief that there is a large excess of 
radiation from these clusters. 

4. The integrated extra-galactic component.—In the brightest part of the 
radio feature the excess radiation (about 80°K at 160 Mc/s) amounts to more 
than 20 per cent of the total intensity in that direction. Most of the background 
emission is, however, galactic in origin. An upper limit on the amount of the 
extra-galactic component has been set by Baldwin (1955) from an analysis of the 
angular distribution of the galactic background. This limit is 500 °K at 81-5 Mc/s. 
Although the spectrum of the extra-galactic component is at present unknown, 
we can use the 81-5 Mc/s observations to show that the excess radiation from the 
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direction of the Supergalaxy is of the same order of magnitude as the total 
extra-galactic radiation. 

If the radiation originates in the Local Supergalaxy, then the first 15 Mpc 
of space measured in a line of sight from the Earth in the direction of the 
Supergalaxy must be intrinsically different from, say, the next 1000 Mpc. 
In view of the large scale uniformity found among the more distant galaxies, 
such a state of affairs seems unlikely. 

5. Conclusions.—The evidence presented in this paper does not support the 
physical association of the radio feature with the Local Supergalaxy and it seems 
possible that it originates within the galaxy. It lies close to the much more 
intense belt of radiation which is very probably galactic. The radiation from 
this latter feature has been discussed by Mrs H. Tunmer (1958). 

Our thanks are due to Mr D. O. Edge for allowing us to use information of 
the 3C survey of radio sources prior to publication. This work was supported 
by a grant from the Department of Scientific and Industrial Research. 
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THE RED SHIFT—MAGNITUDE RELATION IN 
OBSERVATIONAL COSMOLOGY 


W. Davidson 
(Communicated by H. Bondi) 


(Received 1958 August 14) 


Summary 
A reappraisal is made of earlier theoretical work on the red shift—magni- 
tude relation. New formulae are derived which are particularly directed to 
the practical solution of the cosmological problem and the detection of 
evolutionary effects. 





1. Introduction—the cosmological metric.—F ollowing earlier writers we adopt 
as the theoretical basis for cosmological models a space-time description of 
cosmic events characterized by a world metric given by 


dS? = c2 dt®— R(t) 





dr? + r? d@? + r? sin? 6 dd? (a. 

(1 + hkr?/4)? , 
This metric was derived by purely geometrical and kinematic considerations by 
Robertson (1), and also by Walker (2), based on the strict application of the 
cosmological principle. 

The cosmological principle is founded essentially on an assumption derived 
from observations—that as regards the large scale features of the universe there 
is, at all epochs, spherical symmetry spatially about every ‘fundamental ’ observer 
with respect to whom the average motion of matter in his (sufficiently large) 
neighbourhood is zero. In applying the metric (1.1) to the actual universe it 
has to be supposed that local irregularities have been smoothed out and ideal 
spatial isotropy exists about every observer co-moving with the matter in his 
neighbourhood. 

The variable t measures cosmic time as it would be registered on a clock 
carried by such an observer, and R(t) is a disposable function of t whose geometrical 
interpretation is such that k/R*(t) is the uniform curvature of space (¢ = constant) 
atepoch t. ‘The parameter k is a constant taking the values 1, 0 or — 1 according 
as this curvature is positive, zero, or negative respectively. Fundamental 
observers co-moving with the matter of the universe are characterized by having 
constant spatial coordinates r, 0, ¢. 

In observational cosmology the observable quantities employed by astronomers 
are related as far as possible to the theoretical parameters R(t) and k occurring 
in the cosmological metric (1.1). Relations are then developed between the 
observable quantities corresponding to different values of R(t) and k. With 
these the observed relations may be compared and so the means are provided to 
determine those values of R(t) and k which best fit the large scale features of the 
actual universe. 

Since exact relations between the observables do not lend themselves to 
practical tests we shall adopt for the red shift—magnitude relation the method of 
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expansion in series which has been followed by previous writers, notably McVittie 
(3), (12), Heckmann (4) and Robertson (5). For this purpose it is convenient 
to replace coordinate r in (1.1) by means of the substitution 

p=Ry, (1 2) 
where R= R(ty). We take t, to be the present epoch of observation and regard 
it as fixed once and for all so that R, is a constant. The coordinate p will have 
the dimension of distance. The metric (1.1) can now be written 

» R(t) (dp? +p? db* + p*sin® 6 dd? 
dS? = c? dt? — 
0 TD ae 





(1.3) 


where a = 4R,?/k. 
Putting t, —-t=7 we write 
= Ry(t—a,7 + 4ag7?—...) (1.4) 
where ,=R,/Rp, %.=R,/Rp, and so on, a dot indicating differentiation with 
respect to f. 

A light ray in the space-time metric (1.3) is characterized by the equation 
dS? =o, and a photon has the local velocity ¢ relative to a fundamental observer 
of constant spatial coordinates p, 0,¢. Accordingly, the light ray which at epoch t 
is emitted from a source at (p, 8, ¢), and at epoch t, reaches an observer at p=o, 
will by symmetry move along the radial line @= constant, = constant, such that 

t dt a dp 
Oe re il. I. 
BO ole Rew ae 
By expansion of the integrands and integration term by term we find 


er{1 + 40,7 + }(20,?—a,)7?7+...}= 
so that on solving for p we obtain the useful relation 


paer| 1+ dart {i228 24) + Shee...]. (1.6) 


2. Red shift or Doppler effect.—It is now well known that in those models of 
metric (1.3) for which the factor R(t) increases monotonically with epoch a 
shift in wave-length towards the red end of the spectrum will be observed in the 
light from distant sources. This was first shown by Lemaitre (6). In fact, if 
at epoch ¢ a source of constant p, 9, ¢ emits light of frequency v and wave-length A 
according to a fundamental observer at the source, then the frequency v» and 
wave-length A, of the light at epoch ft, when it is registered by another fundamental 
observer, will be given by 


(2.1) 


If the characteristic emissions of an atom are assumed to be of frequency 
independent of epoch, then the observer receiving the light will, if R(t,) > R(t), 
declare a shift of the spectrum towards the longer wave-lengths by a constant 


fraction 
$= a -I, (2.2) 


independent of wave-length. 
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According to the metric (1.3) a progressive increase with epoch of the factor 
R(t) means a progressive increase in the distance between any two neighbouring 
sources of constant p, 6, ¢, if we interpret the metric locally in terms of special 
relativity. ‘The redshifts actually observed by astronomers are now generally 
interpreted in terms of a Doppler effect arising from recession of the sources. 

By expansion of the right-hand side of (2.2) in powers of 7 we can write, 
using (1.4) 

5=a,7+ (4,2 — day)7? + (x3 — xa, + fag)79 +... 


so that, inverting this relationship, we derive 


I hag—a,2\ .. 4 4 24, + da52— da a9\ «. 
r= (—)5+ a-S\s, (5-5 a+3 at — ba1% P4000 GQ) 
Hy a5 a,” 


This result will be employed later. 

3. Astronomical distance, bolometric and heterochromatic magnitudes.—Vhe 
distance normally assigned to a luminous object by astronomers is determined 
by its apparent bolometric luminosity, assuming an inverse square law of diminu- 
tion of intensity and given the absolute power of the source. Thus, if E*(t) is 
the energy radiated at epoch ¢ from the source in unit time over all frequencies, 
while /*)0) is the received energy flux per unit area outside the Earth’s atmosphere 
(before suffering heterochromatic absorption), then the luminosity distance D 
is given by 





E*(t) 
7D (3.1) 


If myo, is the bolometric apparent magnitude of the source, and Mpg) is its 
absolute magnitude corresponding to a luminosity L*po at a distance of 10 parsecs, 
the absolute magnitude of a source being by convention the apparent magnitude 
when at 10 parsecs, we can write 


1*p 0) a 


Myo} “_ Myo = 2°5 logy (L* po} 1* p01). (3.2) 


Hence, if D is in parsecs, 
logy) D=0-2(mbo1 — Moi) + 1- (3.3) 


Since we are unable to measure directly the absolute magnitude Mpo) of a 
very distant source at the time of emission of its light, a relation such as (3.3) 
is not practical for the cosmological problem. However we can observe local 
sources of the same type and obtain at 10 parsecs 


E*(to) 


L* yon = 3° 
47710 


(3.4) 


Therefore, if Mynoy is the corresponding absolute magnitude of such a local 
source, we can write instead of (3.2) 


Myoi — Mynor = 2°5 logio (L*qpon/!*po1), (3.5) 


where myo) and /* 9 refer to the distant source. It follows from (3.1), (3.4) and 
(3.5) that 


bad (4 
log,, D=0-2 {mo — Mynor + 2°5 logy (Fac) +1. (3.6) 
0 
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The introduction of a theoretical luminosity distance in the cosmological 
problem arising from the assumption of an inverse square law of diminution of 
luminosity is due to Whittaker (7). The formula for D appropriate to the general 
metric (1.1) was derived by McCrea (8) who obtained, using the notation of our 
metric (1.3), 


cr Rop 
D= ROG: (3-7) 


Thus D represents the luminosity distance from the origin p=o of a source at 
(p, 4, 6) whose light was emitted at epoch ¢ and is seen at the origin at epoch fy. 

At the time when McCrea obtained this result the depths of space so far 
explored were not sufficient to justify specific concern with the possibility of an 
evolutionary variation of the emission from a typical galaxy. The detection of 
such an effect has now however become of great importance to distinguish an 
evolutionary model (see e.g. Tolman (g)) from the steady state models of Bondi 
and Gold (10) and of Hoyle (11). It is this problem that is our special concern 
in this paper as it seems to the writer that recent treatments are not altogether 
satisfactory since they leave obscurities regarding the practical determination 
of evolutionary effects. 

A formula connecting red shift with bolometric magnitude has been given by 
Robertson (5). This relation is essentially a series expansion of our equation 
(3.6) with D given by (3.7) expanded in powers of 5. It is this formula that is 
adopted by Sandage as the theoretical basis for his analysis of the red shift data 
obtained by Humason, Mayall and Sandage (13). 

The magnitudes directly measured by the astronomers are not of course 
bolometric, corresponding to total energy received outside the Earth’s atmosphere 
(we neglect galactic obscuration as a latitude effect irrelevant to the present 
issue), but ‘heterochromatic’ arising from the selective effects of the atmosphere 
and sensitivity of photographic plates and other detective devices. Accordingly, 
since it has been the custom in the literature to express the theoretical relations 
in terms of bolometric magnitudes, as in Robertson’s formula, a correction has 
had to be applied to the measured magnitudes. 

In the presence of an evolutionary effect, however, this correction as usually 
defined ceases to be directly ascertainable by empirical means, since it depends 
on the evolutionary effects in the spectrum of the source. ‘The process adopted 
by Sandage of allowing for these effects appears to be not well defined. The 
interpretation of the observable data can be made quite clear however by the 
following analysis, by which reference to theoretical bolometric magnitudes is 
made unnecessary. Our treatment is somewhat analogous in principle to one 
given by McVittie (14) but differs in the important matter of dealing with evolution. 

4. Astronomical distance related to photo- or radio-magnitudes.—Suppose that 
from a source whose permanent spatial coordinates are (p, 8, ¢) in the general 
model of metric (1.3) the energy of radiation which is emitted in the waveband 
(A, A+dA) in the interval (t, t+ dt) is 


E(A, t)dddt. 


The function E(A, t) will depend on the type of source considered. This radiation 
will subsequently lie in a spherical shell which will pass over the origin p=o 
in the interval (ft, t)+¢t)), say. At time ¢, the integrated area of the surface: 
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of this sphere according to fundamental observers on it (total proper area) will be 
4np?|(1 + p?/a2)?, 

The total integrated energy of the radiation apparent to such observers will 

only be 
E(A, t) dAdt{ R(t)/ Ry}, 
due to the degradation caused by the recession of the source (eqn. (2.1)). 
Accordingly, the flux of energy per unit proper area at the origin will be 
‘ ( dt t ) 
E(, taal RO 
\dty Ry } 
47p*/ (1 + p*/a?)? 
Since p is constant for the source it follows easily from (1.5) that dt/dt, = R(t)/R,, 
and so the flux of energy becomes 
E(A, t) dX (4:1) 
4nD? ’ 4 





where D is given by (3.7); we may regard D as a useful theoretical parameter. 
To get the energy actually registered on a photographic plate, or other 
detective device, in the received waveband dA, we now require to multiply (4.1) 
by a function s(A,) of recetved wave-length to allow for the selective effect of the 
atmosphere and sensitivity of the plate etc. The energy registered is therefore 


I ; = 
(=) (AE (; ee ) dX 


47D? , 





where we have put A=A)/(1 +4), dA=dA,/(1 +5), in accordance with (2.1) and 
(2.2). ‘Taking account now of radiation emitted over all wave-lengths we find 
for the registered heterochromatic intensity 


I 3 2 
(. = I. S(Ay JE (= ‘ ) dry 


I*net = sab : (4.2) 





We note that if we put s(A,)=1, so that all energy arriving outside the 
atmosphere is registered, we get the bolometric intensity /*po) already given by 
(3.1). ‘This is a consequence of the fact that the total rate of emission E*(t) 
may be written 


E*(t) -| ” E(A, t) dn 


0 


s (5), 8 (2254) (4.3) 


As a particular application of the formula (4.2) we may write down the 
expression for the luminosity of a local source of the same general type, but now 
(conventionally) at a distance of 10 parsecs from the observer at epoch f,: 


| $(Xp)E(Aos to) Ag 
qui0 . (4-4) 





L*ynet) = 
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Here, putting the parameter ¢ equal to ft, allows for the distinction between near 
and distant sources which might arise from an evolutionary effect. Let mnet be 
the registered apparent magnitude of the distant source and Mynet the absolute 
registered magnitude of the nearby source, so that 


mnet — Monet) = 2°5 logig (L* ynet)/!*net). (4.5) 


It follows from (3.5) and (4.5) that the relation between the registered and 
bolometric magnitudes is 

mooi — Maynor = mnet— Manet) — K (4.6) 
where 


(sa) # Ga) 


[7 2004) a 


(=a) [008 (Feat) 


3 $(Xp)E(Apy to) do 


K=2°'5 logio 








—2°5 logio 


=2'5 logi (Eo) 


— 2°5 logis 





(a)fioe(Sy)a] 


[: * sQu)EQo, to) 


noting (4.3). 

The correction to the measured magnitudes which has been made in the past, 
before substitution in theoretical relations involving bolometric magnitudes, is 
the K correction as it is defined in our equation (4.6). It is the correction adopted 
by Sandage, given by his equation B; in the appendix to the paper referred to 
in Section 3. However, Sandage gives an analytic expression for K only for 
the special case of the absence of evolution. This expression, which in our nota- 
tion we shall denote by Ko, is got by putting t= fy in (4.7), so that 


[° OE (85 sto) dy 


. (4.8) 
t. * s(,)EQoy to) Bo 


Ky =2°5 logyy (1 +8) —2°5 logig 





In order to analyse the red shift data by Robertson’s bolometric formula, 
which is equivalent to a series expansion of our (3.6) and takes account of 
evolution, Sandage’s procedure is to make the K, correction, and then to add an 
estimated correction to make up the total K correction appropriate to evolution. 
This procedure therefore does not take account of the fact, as shown by (3.6) 
and (4.7), that part of such an additional correction, viz the term 


25 log 16 laa . 





60 W. Davidson, The red shift—magnitude Vol. 119 


is necessarily cancelled out by the same expression with opposite sign in 
(3.6). This latter term in Robertson’s formula is in fact estimated by Sandage 
separately. ‘The following procedure appears therefore to be more direct and 
employs terms that are all well defined. 

It may be seen that the only relation which it is necessary to consider is that 
got from (3.6), (4.6), and (4.7), namely 


x : A 
| (5) | s(Ao)E (; + 5°) a 
log,, D=0-2| mnet = 


— Manet) + 2°5 logio Px 
| S(Ap )E(Ao, to) 4g 


“ 





TI, 


(4-9) 
in which no reference to bolometric magnitudes appears. Employing the 
relation (4.8) for Ky, which may be regarded as the orthodox magnitude correction, 
this being ascertainable from sources of the required type for which 60, we 
further simplify (4.9) to the form 


| y s(Ap)E (254) dd, 


X 
s(Ag)E (5 ‘) dry 


login D=0:2} mnet — Monet) — Ky + 2°5 logis 





+ ft. 


? ax 


0 


J 

(4-10) 

‘The first three terms on the right are directly derived from observation and 
empirical calculation. ‘The fourth term is an unknown quantity which would 
vanish in the absence of evolutionary trends in the observed spectra. It may 
be expanded, as shown in Section 5, in powers of the red shift 5. The expression 
for D, given by (3.7), will also be expanded in powers of 5 with coefficients 
depending on the parameters «,, %,... of the cosmological model (1.3). The 
equation (4.10) therefore constitutes a direct theoretical relation between two 
observable quantities—registered apparent magnitude and red shift. Comparison 
with observation will therefore provide immediate relations between the unknown 
parameters of the best fitting cosmological model, and also information regarding 
its evolution, if any. 

The result corresponding to our equation (4.10) which is obtained by 
McVittie in his treatment of the problem mentioned in Section 3 may be regarded 
as less satisfactory for the following reasons. In his formula (14, eqn. 8.614, 
p. 157) McVittie separates the total ‘correction’ term into two parts, one that he 
regards as essentially a red shift correction since it does not vanish even in the 
absence of evolution, and the other depending on evolution which is termed 
the Stebbins-Whitford correction. This latter term, however, is misleading, 
since this part of his correction could be non zero even if the Stebbins- Whitford 
effect is zero. ‘This is because the observations of Stebbins and Whitford 
would detect only an evolutionary change in the relative intensities at different 
wave-lengths and not an evolutionary change in the absolute intensities. More 
important is the criticism that McVittie’s redshift correction contains his 
evolutionary parameter at the value corresponding to the epoch of emission, and 
is therefore not a term that can be evaluated directly from observations of local 
spectra, for which we can put ff). This term is in fact expanded by McVittie 
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in powers of 5 whose coefficients correspond to the epoch of emission and there- 
fore, in virtue of our equation (2.3), must themselves be regarded as unknown 
functions of 6. 

The advantages of separating out everything that can be evaluated empirically, 
as in our own analysis, and leaving terms that are non-zero only in the presence 
of evolutionary effects, need to be emphasized. 

5. The red shift—magnitude relation in series form.—The theoretical expression 
for the quantity D which occurs in our formula (4.10) is given by (3.7). Accord- 
ingly, using (1.4) and (1.6), D may be expanded in powers of 7, and in turn in 
terms of 5 by (2.3). Thus, after some calculation, we derive for logy) D in 
terms of 5, as far as the term in 5, the expression 


ras) 1-086 7 _5% ,3/%)\* 
login D = log (< =) + | (+5 23) 3+ {- = 6 a,2 + 4 Na? 
pee > S* I, I 
3 3} ] (5-1) 


It is to be noticed that the coefficient of 5? involves the third derivative of the 
metric factor R(t) for the epoch ty, as well as the curvature of space at this epoch. 
Thus any attempt at finding the curvature from the red shift—magnitude relation 
alone is not to be contemplated. Indeed, in view of the present degree of 
uncertainty in observations, any analysis of the individual contributions to a 
term in 5? is likely to be of little value. Nevertheless if, as is likely, the range of 
5 to which red shifts are measured can be extended beyond the present limit of 
50-2 a term in 5 might well be formally allowed for in fitting the observational 
data to the theoretical relation. This would increase the accuracy obtained for 
the important coefficient of 5. Having remarked this we shall now carry our 
theoretical expansions to order 6 only. 

For the fourth term on the right hand side of (4.10), which depends on the 
epoch t of emission unless evolutionary effects are absent, we write 


51 he nO) (25 2 5") watts © * sQo)é (—% “e sta) B+ O(2) 
25 logyo | é s(Ap)E (25.4) dr 


| * s(Qo)E ota) Do 
| * sQo)E Qos) Bo 








| Beg 2 


+ O(7?) 


| fs yBn ta) dy 
| a1 |” sQo)EQon to) 





8+0(8%) |. (5.2) 





In this we have used the expansion of 7 in terms of 5 given by (2.3). The 
notation E(Ay, ty) means 


75, (EC en eee 
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A point of procedure must now be raised regarding the Ky term in (4.10). 
Since Ky can be evaluated completely by examination of typical local spectra 
we might forego the expansion of this term in powers of 6 and apply the complete 
correction to the apparent magnitudes before fitting the data to a linearised 
red shift—magnitude relation. This was in fact the procedure adopted by 
Sandage in his analysis of the data of Humason, Mayall and Sandage. Alter- 
natively, we may include in the theoretical red shift—magnitude relation the 
linear part of the K, term only and apply the uncorrected data directly to this 
relation, allowing for the coefficient of the linear part of Ky afterwards as a known 
quantity. If a term in 5* were formally allowed for in fitting the data the 
coefficient of 5* would of course be different in the two methods. We shall 
derive the two formulae for these two methods. 

Method (i). In this method the whole of the A, correction is made first to 
each magnitude so that by (4.10), (5.1) and (5.2) the appropriate red shift— 
magnitude relation to order 6 is 


Mnet — Ky= 5 logy (cd) + 1°086 (: + + +s) 8+ Monet)— 5 logy%,—5, (5-3) 
1 


where 
| ‘ S(Ay)E (Ag, to) dX. 
Ae 2 


) 


(5-4) 





nd 


Oy | 5(Ag)E(Ao; to) do 
Method (ii). For this method we require the linear part of Ky. It is easily 
derived from (4.8) that this is 
© } 
[sod Qon to) +208’ Qo to} 
1086 { — - 5, 
[7 :Q0)B0 to) dX 





where 
é 
E’ (Ags bo) = 5 {E(or by )heg const 
( { ) 0Ay ( ' o) ? c 


Accordingly, the appropriate red shift—magnitude relation is now 


Mnet = 5 logy, (cd) + 1-086 (: + — +K* +A) 5+ Monet)— 5 logig%,—5, (5-6) 
1 


where 


| * s(Ap){E(Aos fo) + AoE” (Ags fo)} dy 
Ke= = - , (5.7) 
| s(Ay)E (Nos to) dy 


) 





6. Stebbins-Whitford spectral analysis.—The «* term in (5.6) arises directly 
from the selective effect of red shift and would occur even in the absence of an 
evolutionary process. ‘The A* term would be non zero only in the case of an 
evolutionary universe. It should be noticed that if s(A)) is such that only radiation 
in a narrow waveband is registered on the detective device, e.g. by the use of 
colour filters, so that s(Aj)=o outside this waveband, then the expressions for 
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«* and A* can be written, with an approximation depending on the width of 
waveband employed, 


* ={E(Ao, to) + AgE’ (Ags to)}/ Eo» to), (6.1) 
A* = E(Ag, to)/{ayE (Ao to)}; (6.2) 
for the neighbourhood of any given wave-length Ay. It follows from (6.1) that 


if Ay is chosen so that 
E(Ags to) + AoE’ (Ap, to) =, 


then K*=0, 


(6.3) 


It is clear that this latter condition is likely to be most nearly satisfied at the red 
end of the spectrum, since the peak of the light curve occurs at shorter wave- 
lengths so that E’(Ao, ty) <0 in this region. 

The special form of the expressions for «* and A* given by (6.1) and (6.2) 
would be particularly applicable to the multi-colour measurements of the type 
performed by Stebbins and Whitford (15, 16). It might also find useful appli- 
cation to reception at radio wave-lengths which is usually in narrow wavebands 
(17). 

To deal with the important Stebbins-Whitford type of spectral analysis by 
means of our formula (5.6), suppose that by means of filtering the radiation from 
a distant source is confined to a narrow waveband for which the effective wave- 
length of received energy is Ay. Let the corresponding registered magnitude 
(heterochromatic) be m, and suppose that the registered absolute magnitude for 
the same waveband of a nearby source (60) of the same type is My. Then, to 


order 5, equation (5.6) applies with «* and A* given for narrow wavebands by 
(6.1) and (6.2). Let all these quantities, with dashes, refer to another narrow 
waveband elsewhere in the received spectrum. ‘Then by subtracting the two 
equations written in accordance with (5.6) we obtain 


m—m! = My— My’ + 1-086{(«*—«*’) + (A*—A*’)}S. (6.4) 


Thus, if C is the colour index for these two wavebands in the case of the distant 
source and C, that for the local source, then 


C= Cy + 1°086{(«* — «*’) + (A*—A¥)}5. (6.5) 


The effect of the red shift on the colour index C can easily be calculated from 
a knowledge of «* and «x*’ derived from the specimen source near the observer. 
An additional effect on the colour index, arising from a differential effect with 
respect to wave-length of a systematic evolutionary change in the spectrum, appears 
ifA*—A*’40. Equation (6.5) determines A*—A*’ for different pairs of wavebands. 
This formula appears to be a new one for Stebbins-Whitford type analysis. 
It does not depend on the assumption that the source radiates as a black body. 

By seven colour photoelectric analysis (16) Whitford now finds no significant 
difference in the spectra of elliptic galaxies at varying distances. However, as 
was remarked in Section 4, this type of analysis provides only relative measures 
of the energy in different colours. Until we possess a knowledge of the absolute 
value of A* for one colour in such spectra we cannot say that the total energy is 
not changing with epoch. Nor indeed, until we have such knowledge, can we 
accurately deduce a value for the cosmological parameter «, from the red shift— 
magnitude relation (5.3), or (5.6). As will be shown in a later paper the value 
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of A* can however be derived from the number count data, or a combination of 
number counts and red shift—magnitude data. 

7. Determination of the cosmological parameters «, and %,.—The observational 
data in the form of registered (heterochromatic) apparent magnitudes and 
corresponding red shift can be fitted to our formula (5.3), or (5.6), and the constant 
coefficients found by least squares. As remarked earlier it might lead to greater 
accuracy if a term in 6? was formally allowed for in this process, even if we cannot 
hope to analyse this term profitably. 

From the relation (5.3), following Method (i) of Section 5, we should therefore 
derive values for the constants 

% 
1086 (1+ 34 +2), (7.1) 
and Monet) — 5 10819 % — 5; (7.2) 
while from (5.6), following Method (ii), we find the constant 


1-086 (: + 4 +K* +s) (7-3) 


1 
together with the value for (7.2). ‘The value of «* would be given by computation 
of the expression (5.7) from detailed analysis of typical local spectra and a know- 
ledge of the absorption and sensitivity function s(Ay). We deal effectively, 
therefore, with formulae (7.1) and (7.2). 

Given the registered absolute magnitude Mynet) in the appropriate range of 
wave-length, of local specimen sources of the general type being observed, we 
can arrive from (7.2) at the value of the parameter «,. Since the first approxi- 
mation to the luminosity distance—red shift relation (5.1) is seen to be cd =«,D, 
we observe that «, is the Hubble constant of the expanding universe. The value 
of x, is very sensitive to the value assigned to Mynet). 

Substitution of the value derived for «, into (7.1) would lead to the value of 
the parameter a, if we knew the value of the evolutionary parameter A* for the 
present epoch ty. By identifying the cosmological metric (1.3), in the neighbour- 
hood of any fundamental observer, with the metric of specal relativity, we see 
that «,(=R,/R,) is proportional to the rate at which two cosmological sources 
(o=constant) in the neighbourhood are accelerating away from one another. 
‘Thus the sign of «, determines whether the expansion is speeding up or slowing 
down at the present epoch. 

The expression for A* given by (5.4) can be put in terms of the rate of change 
of absolute (heterochromatic) magnitude of a typical local source at the present 
epoch, thus: 


. 1 f Po ; ) 
Mynet) = — = [ 2s logo 1 | 5(Ag )E (Ao; to) a} | 
0 ( 
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2 Sea ) ] 
=—2°5 = Of10€ 


| : s(Xo)E(Agy to) dg 


= —2'5a,A* logige. 


| * s(Xo)E(Ags to) dp 





Whence At = — 0-92 Mynet)/}- (7-4) 
As mentioned in Section 6 a value for A* could be obtained with the aid of 
the number counts. Meanwhile, the sign of «, for an evolutionary universe 
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could perhaps be deduced if we were willing to assign limits of plausibility to 
Monet. ‘This course is complicated, however, by the possibility of systematic 
selection error. ‘This important matter will be dealt with in the next section. 

8. Observational selection effects.—In the foregoing analysis of this paper the 


possible dependence of the total power of a source, viz., J a E(A, t) dA= E*(t), 


on the epoch t of emission has been attributed to a possibly evolutionary character 
possessed by the source. However, the equations (5.3), (5.6), could equally well 
apply to an apparent evolutionary effect such as would arise by a selection bias 
by the observer which was a function of the distance of the sources and in favour 
of progressively more powerful sources. This is a matter of crucial importance 
for the assessment of a genuine overall evolution in the universe. 

The most obvious case of selection bias that might arise is as follows. If the 
absolute luminosity distribution of the sources is a wide one for every epoch ?, 
then the greater the distance of the detected source the more likely it is that it 
belongs to the brighter end of the distribution, because of the limitations of the 
instruments. There are two cases to be considered. 

If the universe is not in a steady state and there is no creation of matter then, 
in accordance with the usual interpretation, all galaxies were created at approxi- 
mately the same epoch. If with increasing distance the observer selects sources 
intrinsically more powerful than the standard source near him, when at the same 
epoch of emission, then this would violate the assumption that we are following 
the (backward) evolution of the standard source. Since the standard source 
provides the constant Mynet) in (5.3), or (5.6), it is clear that this selection effect 
would have the same influence, via the parameter A*, on the coefficient of 5 found 
from the data as a genuine evolution of the standard source. Consequently, 
if \* were merely estimated by ascribing plausible limits to evolution, instead of 
being deduced independently with the aid of the number counts etc., it would 
lead to a value for «,/x,? lower than the true value. This pseudo-evolutionary 
effect might or might not be comparable with, and in the same sense as, that due 
to a genuine evolution in the time lapse 7( = ¢, — t) between emission and reception 
of the radiation, but in any case the real situation would be obscured. 

If on the other hand the universe is in a steady state then there are two possible 
selection effects. ‘There is the effect described above arising from a possibly 
wide scatter in absolute magnitude of the sources as they are continuously formed. 
In addition there is the possibility that once formed the sources increase in 
brightness as they grow older. In this case the observer would select the oldest 
galaxies at great distances as these would be the brightest. This could lead to a 
very large apparent evolutionary trend in the universe as a whole. 

Hoyle (11) has pointed out that the oldest galaxies in the observable region 
of a steady-state universe will have an age between 10" and 10" years. In fact it 
can be shown that the oldest galaxies indicating a red shift 5 will have an age 


cT8 
Plog. laren 


where 7( = 1/a,) is the reciprocal of the Hubble constant and d is the average dis- 
tance between galaxies where 6 is small. Taking d= 10’ It. yrs. and T = 8 x 10° 
yrs. we find, even for 5 = 0-2, the greatest red shift so far measured, this age is 
4x10yrs. On the other hand the lapse of time, 7, between emission and 
reception is Tlog, (1 +5) which at 6 =0-2 is only 1-5 x 10®yrs, Thus it is possible 


5 
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that in a steady-state universe an apparent general evolution might exist much 
greater than in any genuinely evolving universe. Accordingly, if this were not 
allowed for by an independent determination of A*, a value for «,/«,” too small, or 
even negative, would be deduced so that the steady state model, on this test at least, 
would be inadvertently ruled out. 

In the case of the first type of apparent evolution, arising from a scatter in 
absolute magnitude, a basis for the estimation of this effect in both kinds of universe 
might be found from a study of the luminosity function of galaxies showing small 6. 

We point out that the presence of either type of selection effect could be detected 
by the use of measurements of angular distance between distant galaxies in con- 
junction with red shift. It will be shown in a later paper that there is a simple 
relation between angular distance and red shift for two ‘average’ galaxies showing 
this red shift, either in a steady-state universe or in a genuinely evolving universe. 
For instance, in the case of the steady-state universe, if d is the average distance 
between galaxies where 6 = 0 (i.e. where the galaxies are relatively young), then 
the angular distance, 6, between two galaxies which are at the same distance d 
apart, but both showing sensible red shift 5, will be given by 


p= 4 (+8) 
crt 6 


Consequently, if it is found that the galaxies observed are on the average at a 
greater angular distance apart, this would mean that the observer is viewing 
galaxies either older than the average age of the galaxies in his own neighbourhood 
or at least rarer, and therefore probably brighter, than the galaxies of his 
neighbourhood. 

The essential features of this paper formed a chapter of a Ph.D. thesis 
submitted to the University of Londonin March, 1958. The author is grateful for 
permission to publish it. 


Royal Holloway College, 
Englefield Green, 
Surrey: 
1958 August 12. 
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A PROPOSED ASTRONOMICAL TEST OF THE “BALLISTIC”’ 
THEORY OF LIGHT EMISSION 


Herbert Dingle 
(Received 1958 October 24) 


Summary 


Ritz’s ‘‘ ballistic ’’ theory of light emission is outlined, and the evidence 
on which it has been regarded as untenable is re-examined. It is shown 
that the theory as presented by Ritz is indefinite, and that, when it is amended 
to a form consistent with the postulate of relativity of motion, it is no longer 
vulnerable to the earlier criticisms. "The question of its admissibility is 
therefore open. It is pointed out that determinations of the ‘‘ constant 
of aberration ’’ from observations of distant nebulae might settle this ques- 
tion, and such determinations are recommended. 





It is well known that, according to Einstein’s restricted theory of relativity, 
no material body can move with a greater speed than that of light in vacuo, 
namely, about 3-10!° cm/sec. Recent work in nuclear physics, however, seems 
to conflict with this*, for, according to the uncertainty principle, the uncertainty 
in the velocity of a body may approximate to infinity. Independently of this, 
the internal consistency of the restricted relativity theory seems questionable 
if the postulate of the constancy of the velocity of light is given its usual inter- 
pretation of requiring that all bodies present at the origin of a wave of light, 
whatever their relative motions, must be regarded as remaining thereafter at the 
centre of the light-spheret. These difficulties are removed if the postulate 
be interpreted merely. as requiring that the velocity of light relative to its actual 
material source shall always be c, and this hypothesis was, in fact, advanced 
long ago by Ritz{ and was strongly supported by la Rosa among others: it is 
usually known as the “‘ ballistic’ hypothesis of light emission. ‘Though for a 
time regarded as a serious rival to Einstein’s theory, it was ultimately rejected as 
being contrary to observation, and has now for many years been held to have 
been definitely disproved. In view of the new difficulties facing Einstein’s 
theory, however, a revision of the hypothesis and of the grounds for its rejection 
seems desirable. 

Ritz’s hypothesis is based on the postulate of relativity of motion, which 
he expresses in the following way§:— “‘ La seule conclusion qui... me semble 
possible, c’est que /’éther n’existe pas, ou plus exactement, qu’il faut renoncer 
a se servir de cette image; que le mouvement de la lumiére est un mouvement relatif 
comme tous les autres, que les vitesses relatives seules jouent un réle dans les lois de 
la nature.” 

On this basis he makes the following statements concerning the motion of 
light\|: “‘ les particules lumineuses expulsées en tous sens a l’instant t se meuvent 
avec une vitesse radiale constante et remplissent constamment une sphére dont le 

* See, for example, W. Heisenberg, The Physicist’s Conception of Nature, p. 48, London, 1958 

t Bull. Inst. Phys. 9, 314, 1958. 
t+ W. Ritz, Ann. de Chim. et de Phys., 13, 145-275, 1908. 


§ Ibid., p. 207. 
|| Zdid., pp. 210-1. 
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centre est animé du mouvement de translation w qu avait P a l instant de Il’ émission; 
st w est constant, ce centre continuera donc de coincider avec P.... La vitesse 
de la lumiére dépend donc de celle que posséde le corps qui l’émet au moment 
de l’émission; a partir de cet instant, la vitesse des particules reste invariable, 
quel que soit le mouvement ultérieur de P.” 

(It should be remarked that, although Ritz speaks of “light particles ”’, 
his theory on this point is purely kinematical and implies nothing whatever 
concerning the nature of light. As Whittaker has pointed out*, a decision on 
the “ballistic” hypothesis is “ without significance one way or the other in 
the dispute between the wave and corpuscular theories of light ’’.) 

Unfortunately, this statement of the hypothesis is indefinite, since Ritz 
fails to make clear to what standard of rest the velocities are referred. ‘The 
ether having been abjured, the phrase, ‘‘ mouvement de translation w qu avait 
P a Vinstant de lémission”’, leaves us guessing with respect to what P (the 
source of light) has the velocity w; and the phrase, “la vitesse des particules 
reste invariable, quel que soit le mouvement ultérieur de P’’, seems to indicate 
that, whatever the standard may be, P may have a varying velocity with respect 
to it, and therefore with respect to the light particles after the moment of their 
emission. De Sitter, who has given what is generally considered the most 
convincing disproof of the “ ballistic ” hypothesist, tacitly assumes the standard 
of rest to be the Earth (or perhaps the Sun). He considers observations of a 
binary star—for simplicity let us suppose the plane of the orbit to pass through 
the Earth and one component only to be visible. If the orbital velocity is z, 
the light issuing while the star is approaching the Earth will have velocity c+, 
and that issuing while it is receding from the Earth will have velocity c—v. 
If D is the distance of the star from the Earth, the former light will take time 
D (c+v), and the latter time D (c—v), to reach the Earth. The orbit calculated 
from the observations will therefore be very different from the actual orbit, 
and if the latter is Keplerian the former will not be so, but, as de Sitter showed, 
may in a particular case imply that the star has three different velocities at the 
same time. ‘The actual observations, however, do correspond to a Keplerian 
orbit. We must therefore, said de Sitter, reject the supposition that the velocity 
of the light depends on that of its source. 

If we are entitled to suppose that a member of the Solar System is acceptable 
as a unique standard of rest, then de Sitter’s evidence against the “ ballistic ”’ 
hypothesis is final. But a different interpretation of Ritz’s theory was given by 
Tolman{, according to whom the theory requires that ‘‘ throughout its whole 
path light retains the component of velocity v which it obtained from the original 
moving source. ‘Thus all the phenomena of optics would occur as though light 
were propagated by an ether which is stationary with respect to the original 
source.” Qn this view, therefore, the light particles (or wave-front), at a time t 
after emission, would be at a distance ct from the source, no matter how the 
source might have moved, between emission and the instant ¢, with respect to 
any standard at all. 

Ritz regrettably died shortly after his paper appeared, and he therefore had 
no opportunity of elucidating his meaning. It is clear, however, that the 

* E. T. Whittaker, History of the Theories of Aether and Electricity, Vol. 11, p. 39, London, 1953. 


t W. de Sitter, Proc. Amsterdam Acad., 15, 1297, 1913; B.A.N., No. 64, 2, 163, 1924. 
¢R. C. Tolman, Phys. Rev., 35, 136, 1912. 
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only “ ballistic ” hypothesis that is consistent with the relativity of motion, the 
absence of a luminiferous ether, and any plausible form of postulate of constant 
light velocity, is that contemplated by Tolman. Suppose that a body, alone 
in space, emits a pulse of light. If there is any meaning at all in calling the 
velocity of light constant, the pulse must continue to move at velocity ¢ with 
respect to the source. In the absence of a universal stationary medium, there is 
nothing else with respect to which to express its velocity. It will make no 
difference to the relative velocity of the light and its source if we place a body some 
light-years away with respect to which the source moves to and fro; so if, from 
the point of view of that body, the velocity of the source fluctuates between +v 
and —v, the velocity of its light must fluctuate between c+v and c—v, and a 
beam issuing ¢ seconds after an earlier one will reach a distance D from the 
source precisely ¢ seconds after that earlier one. 

If, then, we interpret the “ ballistic’? hypothesis in this way, clearly the 
observations cited by de Sitter afford no evidence against it. With respect to 
the Earth, each particular light-pulse issuing from the star fluctuates in velocity 
during its journey, in unison with the fluctuations in velocity of the star, and 
the observations therefore give precisely the same relation of orbital velocity 
to time as that characterizing the actual orbit. Such fluctuations seem strange 
only if the idea of a universal ether remains unrecognized in our minds. If we 
consider only the Earth and the star, the apparent movement of the star will be 
indistinguishable in character from that which we attribute to the effect of 
aberration of light. It is only when we take the rest of the universe into account 
that we realize the convenience of ascribing the first to the orbital motion of 
the star and the second to the orbital motion of the Earth. It is tempting, of 
course, to think that the overwhelming measure of this convenience must make 
the whole material system of the universe a preferred standard of rest, but 
however that may be in general mechanics (it is unnecessary to consider that 
question here), it is an impossible supposition in this case, for we are concerned 
only with the velocity of light in space. If we suppose that the “ universe ” 
provides an absolute reference frame for that, then the Michelson—Morley 
experiment should have shown the Earth’s velocity with respect to this frame, 
but it did not. Accordingly, ‘‘ universe’ becomes in this connection merely 
another name for “ luminiferous ether”, and the inapprehensibility of that 
leaves us with the source of light as the only possible reference body with respect 
to which the velocity of light can be postulated to be constant. 

It should be unnecessary to add that, on the “ ballistic’ theory, the simple 
addition law of composition of velocities holds, so that the velocity of light with 
respect to the receiving body may exceed c. 

Tolman, however, went on to propose a different test of Ritz’s hypothesis. 
He held that the hypothesis would require a positive result in the Michelson- 
Morley experiment if the light used came from an astronomical source moving 
with respect to the Earth, for in that case the velocities of the light along the 
two arms of the apparatus would be different, and a change of fringe system 
would be expected on rotating the apparatus. Later, the experiment was 
repeated with celestial light by Tomaschek*, and independently by Millert, 


* R. Tomaschek, Ann. d. Phys., 73, 105, 1924. 
+ D. C. Miller, Proc. Nat. Acad. Sci., 11, 306, 1925. 
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but with the same result as in the original performance. An experiment by 
Majorana* with moving terrestrial light also gave a null result. ‘Tolman accord- 
ingly concluded that the “ ballistic”” hypothesis was disproved. 

It appears, however, on general grounds that experiments of this kind can 
afford no evidence on the matter, for, since the only velocity concerned is the 
relative velocity of the apparatus and the incident light, an experiment which 
failed to reveal it for one value could hardly do so for another. The equality 
of the two paths taken by the light in the Michelson—Morley experiment is 
independent of the velocity of the light when it enters the apparatus. If the 
velocity should change, the time taken by the light must change equally along 
both paths, and no differential effect would be expected to be observed. In 
terms of interference fringes, it should be noted that on the “ ballistic’’ theory 
the frequency of the (monochromatic) light received is proportional to the velocity 
of its reception. A change in that velocity would therefore show itself as a 
change in frequency, and in no other way. But a similar change of frequency 
would occur with a change in velocity of the radiating body on any theory: it 
would be observed as a Doppler effect. Consequently, any change in the 
fringe system arising from a change in the velocity of the radiating body relative 
to the apparatus would not distinguish between the “ ballistic’’ theory and 
Einstein’s theory of invariability of velocity of light reception. 

In view of these general considerations, it is necessary to examine 'Tolman’s 
argument to see where the discrepancy lies. He supposes the experiment 
performed with light from the Sun when that body is moving along the direction 
of the arm AB with velocity v with respect to the apparatus, and he writes:— 
“It is easy to see that the Ritz theory would lead us to expect c+w for the 
velocity of light in the direction AB, c—v for the velocity in the opposite direc- 
tion.” But the theory leads us to expect nothing of the kind. The velocity 
of the light in the direction AB is certainly supposed to be c+v with respect 
to the apparatus, but the velocity in the opposite direction depends on the law of 
reflection of light. On this the Ritz theory says nothing, and the only reasonable 
assumption is that on this, as on all non-ether theories, light is reflected as a 
perfectly elastic body, and leaves the mirror at the same velocity with respect 
to it as that which it had on incidence (but of course with opposite sign). In 
that case the shift of fringe system assumed by Tolman would not be expected 
to occur. 

This appears to exhaust the evidence on which the “ ballistic” theory has 
been rejected, and we must conclude from these considerations that at present 
there is no basis at all on which to decide for or against the theory as interpreted 
by Tolman, namely, that the velocity of light proceeding from a material source 
remains constant, at a value c, with respect to that source, and is received at 
velocity ¢+v (vectorial addition) by a body moving at velocity v with respect 
to the source at the time of reception. The purpose of this note is to suggest 
how legitimate evidence might be obtained. 


The “ constant of aberration” is, to a first approximation, proportional to 
the ratio, v/V, of the orbital velocity of the Earth to the velocity of the light 
received from the celestial body concerned, both measured with respect to the 
Sun as a standard of rest. On the “ ballistic”? theory this quantity should 


* Q. Majorana, Phil. Mag., 37, 145, 1919. 





No. 1, 1959 test of the ‘‘ ballistic’’ theory of light emission 71 


vary with the radial velocity of the celestial source, while on the view that light 
is always received at velocity c, it should be constant. Stellar velocities are too 
small to distinguish between these possibilities, but a distant nebula might move 
fast enough. ‘The light from a nebula with a velocity of recession of 30 000 km/sec, 
for instance, should, on the “ ballistic”? theory, be received at a velocity 
o-gc, and the aberrational constant should accordingly be about 23” instead of 
20"-5—an easily detectable difference. Still larger differences are within range 
of detection. 

Unfortunately, if such a nebula yielded the normal value for the constant 
the result would not necessarily be conclusive, owing to the uncertainty which 
still remains as to whether the nebular red-shift is properly attributed to reces- 
sional movement. On the other hand, if the larger value were obtained it would 
be difficult to think of any possible explanation except that the red-shift is due 
to recession and the “ ballistic’? theory is correct. Furthermore, a method 
would become available of finding the recessional velocity of any nebula bright 
enough to be directly photographed. These possibilities make it very desirable 
that this observation should be made. 


104 Downs Court Road, 
Purley, Surrey: 
1958 October 23. 





NOTE ON 'THE MEASUREMENTS OF RADIO EMISSION 
FROM CASSIOPEIA A 


V. Pletchkov and V. Razin 


Measurements of radio emission from the discrete source Cassiopeia A at 
a frequency of 300 Mc's have been attributed to us by a number of authors 
(1, 2,3). We wish to point out that this is erroneous. We have made the 
measurements of intensities of radio emission from the Cassiopeia A and two 
more discrete sources at wave-lengths of 9-7 cm (frequency is 3100 Mc/s) (4) 
and 3-2cm (4,5) only. It should be noted that the flux density of radio 
emission from Cassiopeia A cited in the first three references for a frequency 


of 300 Mc s is very close to that measured by one of us (Razin) at a frequency 
of 3100 Mc's, viz. 1-3 x 10-73 w/m? (c/s) (4). 
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PESIOME JIOKJIAJIOB, B TTEPEBOE HA PYCCKHM A3bIK 


O @PMOPME KOMETHbIX XBOCTOB 
B. JT. Mux 


Ha ocHope mpeycTaBneHHua O KOMeTe, KaK O poe YaCcTHI, OOpasoBaBUIMXCH BHYTPH KOMCTHbI 
BONM3SM NepHresIMA, BbIUHMCIeCHO MpOCTpaHCTBeHHOe paciipewesyeHue waCTHI, pasHbIX pasMepoB 
UIA prlia NOSHeEHUIMX MOMCHTOB. Ornenenue HCKOTOPbIX 4YaCTHI OT rlaBHOM 4wacTH KOMeTbI 
OOyCNOBNeHO CBeETORbIM aBJICHHeM (HMJIH JMHaMHUeCCKH SKBHBaJICHTHbIMM 3¢pdeKTaMM) MH MO?KeT 
ObITh TpefcraBlreHO Kak ocnaGmeHHe COJHeCYHOTO NPHTsPKeHHA, BKIIOUaA, KOHeUHO, jaKe 
riepeMeny 3Haka. OOcy2KaloTCA HEKOTOPbIe CBOMCTBa TeOMeTPHU4eECKOrO MeCTa YACTHI (C pasHbIMH 
pasMepaMH HM pa3SHbIM 93¢p(PeKTHBHbIM MIPHTAXKeHHeM), OJHOBPeCMeHHO MOKHalOuMx KOMeTY; 
HeKOTOPble H3 MOJYYCHHbIX Pe3yIbTAaTOB MPHMeCHAIOTCA K YaCTHOMY CJIy4alO KOMeTbI ApeHja- 
Ponana. Teopna yaetT mpamoe oOObACHeHMe BHMMOrO OCTporo BBICTyMa, HallpaBeHHOrTO K 
Connuy, KoTopbiii HaGmojaicA y STOM KOMeTbI, KaK Ppe3yIbTaTa BeepooOpasHoro pacnpeyleseHHA 
OT CIMBINHXCA UACTHLL, KOTOPOe MO?KET HaGIOMAaTbCA TOIbKO CCIM PAaCCMATPHBATb CTO M3 NYHKTOB 
nevKall{Mux B ero MIOCKOCTH HIM O4eHb 6M3KO K Hell. 


YVINHPEHHE BCJIEJKCTBHE CTOJIKHOBEHHH HU C]IBHT 
PESOHAHCHONM JIMHHH KAJIBILNA 


B. P. Xatindmapu 


Hamepenbi CIBHT BCIeXCTBHe CTONKHOBeHHH HM yuIMpeHHe IMHHii Ca 4227 A, o6ycnoBneHHBIe 
BHCIUHHM jlaBlieHHemM renusa. JIMHMA NoNy4alacb B NOPIOWeHHH B pe3syAbTaTe NpoxorxKTeHHi 
Ses10ro CBeTa UEPe3 Napbi KaJIbUMA B MIPHCYTCTBHH PeJIMA NPH PasJIMGHBIX aBIeCHHAX MeHbile I aTM. 
CnektpasbHoe paspemenne Opi0 oGecneueHo pemrerKoii Ba6koka, HCnONb30BaHHOM B HOBOM 
Oxccbopackom cnekTporemmorpade. Tlonyumpwua sHnnn, oOycnoBneHHaA 3aTyxaHHeM, Gbuia 
HaiijeHa paBHoii 1,72+0,12. 10° ** cm | Ha I aTOM Ha I CM* resIMA, a CBHI B CPHONeETOBYIO CTOpOHy 
cocTaBJIAJI 0,0§ +0,04 . 10°*° cm”! Ha I aTOM Ha I cm*. TlooTomy oTHOUIeHHMe yUIMpeHHA K CABHTy 

>19, 4YTO 3HaUMTebHO Gombe BeJINYMHbI 2,76, mpeacKasbiBaemoi Teopneli JImHyzxXombMa. 
Cupur TakoKe, cormacho TeopHn, AOJpKeH ObITh K KpaCHOMY KOHL. 

OGcyKMaIOTCA BOSMO*KHbIeC OOLACHCHHA 9TOTO pacxoneHHA. Tloka3aHo, YTO OTTAIKHBaTe- 
JIHbIe CMIbI OHskoelicrBHA, elicTByIOlHe MexK]{y aTOMaMH Ka@JIbluiMAl M TeIHA JOJDKHBI 
VUHTbIBaTbCA Tak *Ke, KaK HM CHMJIbI JtaibHOselicTBHA BaH Rep Baanpca. M3 usmepennii Ha OCHOBe 
Teopun JIMHaxXobMa OeHeHbI KOod*uUMeHTHI JIA cust O6onx THNOB. B sakmioueHHe MoKa3aHo, 
uTo B ycuoBHAx ComHua KosdpduuMeHT 3aTyxaHHA ya Ca 4227 A menbule 2,2.10-*cm* ana 
+=0,1 H YTO CUBHMr, OOYCIOBIICHHbIM CTOJKHOBCHHAMH, JaeT mpeHeOpexkHMO MaJIbIi BKNay B 
KpacHoe CMeLileHHe COJTHCUHDIX HHH, 


CIIBHTH JIMHH B CIMEKTPE APrOHA BCJIEIICTBHE CTOJIKHOBEHHT 
B. P. Xaiindmapu u K. A. Tomac 


Bbui H3MepeHbI OGYCHIOBICHHbIC CTOJIKHOBCHHAMH CXBHIM {Byx J1MHMii aproHa A 7067 A Hu 
16965 A; AHHHM BOSOy2KTaIMCh NOCTOAHHbIM TOKOM IPH HCKPOBOM pa3pajle B YMCTOM aproHe. 
JIlmHuM M3y4aINCh IPH MOMOLWIH HMHTepdepometpa Dabpu-Ilepo, ckpeuleHHoro Cc NPH3MeHHbIM 
cnekTporpadbom Xwuuprepa. JlwameTpbl MHTepdbepeHUMOHHbIX KOJIeIl H3MCPAJIMCh JIA pAsla 
jlaplennit aprowa MH TOKOB pa3spaya. IlokasaHo, 4TO M3MepeHHbie CBHTH Gblin OOyCOBJICHbI 
CTONKHOBCHHAMH M@KJy ABYMA aTOMaMM aproHa, a He BIIMAHHeM 3aps>KeHHbIX YacTHy. CaBur 
iia 47067 coctaBiseT: —0,61+0,05.10°*°cm™! Ha I aTOM Ha I CM® aproHa, H jUIA A 6965: 
—0,65+0,05 . 10°** cm! Ha I aTOM Ha I Cm’ aproHa; MMHYC COOTBeTCTBYeT CMCLIICHHIO B KpacHbiit 
KOHEIL CUeKTpa. OTH pesyIbTaTbI HAaBOAAT Ha MBICJIb, UTO CHJIbI GAM3KOReCHCTBHMA OTTAIKMBaHHA 
M@?KJY ATOMaMM aproua MOryT GbITb MaJIO CYLLCCTBCHHbI NIpH CABHTax JIMHH, B OTUMUHE OT Ciry4an 
NMHHM Kalba A 4227 A, BOsMyUIeHHOH resem. ‘ 

? 
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HM3MEPEHHA YPIIOBOrTO JTHAMETPA HCTOUHHKOB PATHOM3SIIVYYUEHHA 
JIEBEJ]Ib (19N4A) H KACCHOTTEA (23NsA) HA BOJIHE 10,7 CM 


B. Poycon 


OnwcpiBaeTcsA almapatypa, HCNOb30BaHHaA [UIA H3MePCHHA YIIOBbIX HaMe€TPOB HCTOUHHKOB 
PaqHOHBTy4eHHA Ha BOlIHe 10,7¢M. IIpnHBeyeHbI pesybTaTLI H3MepeHHi IA ABYX B3aHMHO- 
NePNeCHAMKYIAPHbIX HallpaBlieHHit Kak JUIa Kaccnonen (23NS5A), Tak u aia JieGena (19N4A). 
OTH pe3yJIbTaTbI NOKa3bIBalOT, UTO NePpBbIit HCTOUHHK XapakTepH3yeTCAH KpyroBoll CHMMeTPpHeit 
H HMMeeT IIPHMePHO TOT 2Ke PasMep, UTO H Ha METPOBBIX BOJIHAaX, B TO BPeMA KaK BTOPOM HCTOUHHK 
3aMe€THO aCHMMeTPHYeH C GOON OCbIO, KOTOPaH. BO3MO)KHO, HeCKOJbKO Gombe Ha Goree 
BbICOKHX YacToTax. LlosHuMOoHHbI yron Gombulol OCH CocTaBAeT 109° + 2°. 


OJVIMMICOH], CKOPOCTEH H PPABHTALINOHHBIM TIOTEHI[HAJI TAJIAKTHKH 
II. A. Votiman 


B 19s6r. M. II mMuur nocrpons Mojesb PpaBNTalHMOHHOrO NMoTeHUMasa DaswakTHKH Ha OCHOBaHHH 
oO6paSoTKH KPyrOBbIxX CKOpocTeli, NMOJYYCHHbIX H3 21-CM paqHoHaGOWeHH WM Apyrvx WaHHbIx. 
BbIACcHAeTCH, YTO C HeOONbIIHMH H3MeHeHHAMM MoOjIeb IIMuaTa MOXKHO pacCMaTpHBaTb Kak 
NpHHaWierKalllyiO K THITY CTaWUMOHaPHbIX rataKTHK, ONMCaHHOMY B 1941 Tr. Kammom. Osmncony 
ckopocteii B ramaKTHKax 9TOrO THIMa MO7KeT GLITh TPeXOCHbIM H MO?KHO HCCIe{OBaTb H3MCHCHHe 
9THX Oceii B pa3smMUHbIx YacTtAx TanakTukuH. B KauecTBe OAHOrO H3 pesyIbTaTOB mosyuaeTcA 
HaKJIOH ,, JIMHHM BEpTeKCOB ”’ K raslaKTH4eCKO MIOCKOCTH Ha CPaBHHTesIbHO HEOOJbILIOM VAaJIeHHH 
OT 9TOH MIOCKOCTH. 


PAJIMOH3JIYUYEHHE H3 PAHOHA ,, CBEPXTAJIAKTHKH ” 
Jl. P. Metixuwadm u JIac. 9. Boadyun 


HaGmoeHuA MpOBOAWIMCh MPH NOMOUIM KeMOpHDKCKOrO pasMOTesIeCKOMa Ha BOJIHE 1,9 M. 
Ha ocHospe 9THX HaOmoqeHHH ObuIM MocTpoeHbI H30qdorn: obsacTH Heda 10h-14! ¢ wesbIO 
UCCIe€OBaHHA Psa eTaseH, KOTOpbie OLIN NPHIMHCaHb! H3J1yYeCHHHO OT ,, MCCTHOM CBepxTaslakKTH- 
KH”. Koppesiauna pacnpeyeneHHA OMHCCHH C aHHbIMH MoOjCueTa TalaKTHK OKa3aslach He 
yocTaTOuHO TecHoil. HaOmogaemble HHTeHCHBHOCTH B HeCKONbKO pa3 OoNbUIe, YEM MO?KHO 
ObLIO OOKHMLaTh [VIA HOPMAJIbHBIX PaslaKTHK, a MOACUET PasHOHCTOUHHKOB MOKa3aJl, YTO ITOT 
PeSYJIbTaT HEIb3A OOLACHHTh MaJIOH Orel! TanaKTHK, HMCIOWLMX AHOMAJIbHO BbICOKHe CBETHMOCTH 
B paqHoyManasone. PayqMoApKoctTb OHO H3 JeTatei 3HAaYHTeIbHO Gombe, YEM APKOCTb BCCro 
(poHa BHeratakTHYecKorO H3TyueHHA. Bpicka3prBaeTcCA MpeiouoOx*KeHHe, YTO 3TH jleTaIH B 
JlCHCTBUTeENbHOCTH He CBA3AHbI CO ,, CBepxraslaKTHKOH ”’. 


COOTHOIUIEHHE ,, KPACHOE CMEITLEHHE—3BE3]IHAAH BEJIHNUHHA ” 
B HABJHOJIATEJIBHOH KOCMOJIOrTHI 


B. Jletiaudcou 


CyentaHa mepeoueHKa OnyOJIHKOBaHHOHN paHee TeopeTHueCKOH padoTbI MO COOTHOLICHHIO 
» KpacHoe CMellleHHe—3Be3qHaAd BelMuHnHa”. BpiBeseHbl HOBbIe (OpMyYJIbI, KOTOPbie cnemHa- 
IbHO pe {HaSHa4eHbI! UIA MpaKTHYeCKOrO PewICHHA KOCMOJIOrHYeCKOH MpodsieMbI WH OOHapy?*KeHHA 
IBOJHOUMOHHDIX 3¢pcpeKTOB. 


BOSMOXKKHbIN ACTPOHOMHYECKHH KPHTEPHH JIA MPOBEPKH 
» BAJIIMCTHUYUECKOH ” TEOPHH CBETOBOrO H3/IYUEHHA 


X. JTunea 


Kpatko wanaraetca ,, GammcrMueckad””’ TeopuA cBeTOBOrTO H3iyueHHA PutTua HM 3aHOBO 
epeCMaTPHBaloTcH apryMeHTbI, Ha OCHOBE KOTOPbIX 9Ta TeCOPHA CUMTAaaCh HECOCTOATEJIBHON. 
[loka3aHo, 4YTO 9Ta TEOPHA B NpescTaBsIeHHOM Putiuem Bue HeolpewesieHHa, HO eCJIM ee JONOJHUTb 
M MIPHBecTH K (bopMe, COBMeCTHMO)i C MOCTY1aTOM OO OTHOCHTesIBHOCTH JIBH>KCHHA, OHA CTAHOBHTCHA 
He€YASBHMOH {UIA Iperx*KHeli KLPHTHKH. YkKa3biBaeTca, YTO ONpeseeHHA ,, MOCTOAHHOM aGeppaun ” 
H3 HaOmOeHHH VilasieHHbIX TanlakTHK MOrIH Obl pelIHTb 93TOT BOoNpoc. PekomeHjyeTca 
nmpoBesleHhe TaKHX OMmpestesteHHil. 
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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by 2 summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. ‘The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 
ryped using double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 

H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
1.A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I.A.U. practice. 


4. Diagrams.—These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. Tables.—These should be arranged so that they can be printed upright on 
the page. 

6. Proofs.—Costs of alteration exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts.—When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MSS. 


Reading of Papers at Meetings 
8. When submitting papers authors are requested to indicate whether they will be 
willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 
g. Postcards giving the programme of each meeting are issued some before the 
meeting concerned. Fellows wishing to receive such cards whether i 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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